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Abstract. The box-ball system is an integrable cellular automaton on one 
dimensional lattice. It arises from either quantum or classical integrable systems by 
the procedures called crystallization and ultradiscretization, respectively. The double 
origin of the intcgrability has endowed the box-ball system with a variety of aspects 
related to Yang-Baxter integrable models in statistical mechanics, crystal base theory 
in quantum groups, combinatorial Bethe ansatz, geometric crystals, classical theory of 
solitons, tau functions, inverse scattering method, action-angle variables and invariant 
tori in completely integrable systems, spectral curves, tropical geometry and so forth. 
In this review article, we demonstrate these integrable structures of the box-ball system 
and its generalizations based on the developments in the last two decades. 
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1. Introduction 



1.1. The box-ball system 

The box-ball system (BBS for short) is a cellular automaton introduced by Takahashi 
and Satsuma in 1990 [71]. It is a dynamical system of finitely many balls in an infinite 
number of boxes aligned on a line, whose time evolution is given by the following rule. 
We assume that each box can accommodate one ball at most. 

(i) Move the leftmost ball to its nearest right empty box. 

(ii) Move the leftmost ball among the rest to its nearest right empty box. 

(iii) Repeat (ii) until all the balls are moved exactly once. 

This defines an update corresponding to the one time step t — )■ t + 1. We remark that 
the above evolution rule is invertible. Let us show an example. By starting with the 
following configuration at time zero, 



t=0 



we obtain the configuration at t = 1 as 



t=l 



and so on, ... 



t=2 
t=3 
t=4 



One observes that "a series of three balls" and "a series of one ball" proceed to the 
right stably unless they are "too close" to each other. The larger one is faster than the 
smaller one, so they eventually collide. After the collision, it is non-trivial that they 
come back in the very original shape as 3 + 1 — > 1 + 3, instead of being smashed into 
pieces like 3 + 1— J-l + l + l + lor getting glued together like 3 + 1 — )■ 4. Moreover, 
the collision has caused a phase shift; observe that the trajectory of the larger (smaller) 
series has been shifted by 2 to the right (left). 
Let us watch another example: 



t=0 
t=l 
t=2 
t=3 
t=4 
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Here we have three series of 4, 3 and 1 balls from the left, and they are interchanged 
into the reverse order 1, 3 and 4 during 1 < t < 3. These behaviors of series of balls 
remind us of solitons in the theory of nonlinear waves. We call the number of balls in a 
series of balls before or after collisions an amplitude of the soliton. (A precise definition 
of solitons and their amplitude will be given later. See for example f l2.47p or f l3.22p .) 

One can also set up the BBS with the periodic boundary condition [8^. Let L be 
the number of boxes aligned on an oriented circle. We put M < L/2 balls into them. 
The balls are moved by the same rule as the previous (i)-(iii) for the original (infinite) 
BBS except a minor adaptation to the fact that nothing can be "leftmost" on a circle. In 
(i), the procedure can be started from any chosen ball. In (ii), the terms "leftmost" and 
"nearest right" are to be understood along the direction of the orientation of the circle. 
Then the modified evolution rule is well-defined in the sense that the result is actually 
independent of the choice of the first ball to move. Moreover it is again invertible. Let 
us look at an example of L = 12 and M = 4 in the following, where we identify the left 
and right boundaries (thick lines): 



t=0 
t=l 
t=2 
t=3 
t=4 
t=5 
t=6 



] 



One can observe that the larger soliton overtakes the smaller one repeatedly. The biggest 
difference from the infinite BBS is that the system has now a finite configuration space; 
there are just (^) states. Thus, any state is cyclic, i.e. by starting with an arbitrary 
initial state one comes back to itself in a finite time. 

Let us motivate our study on BBS from the viewpoint of solitons and integrability. 
Nowadays, the term soliton is widely used to mean, somewhat loosely, various special 
solutions to nonlinear equations that exhibit particle like behavior or certain stability. 
In its original context of Zabusky and Kruskal [85] however, it meant a solitary wave 
solution in an infinite dimensional nonlinear dynamical system (KdV eq. mentioned 
below) with more stringent properties as follows: 

(a) particle-like propagation (constant velocity, stability under multi-hodj collision), 

(b) factorization of scattering (pairwise scattering with phase shifts). 

Existence of solitons is a signal of integrability, which at least postulates an infinite 
number of conserved quantities (integrals of motion). The historically important and 
famous integrable systems of such kind are Korteweg-de Vries (KdV) equation and 
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Kadomtsev-Petviashvili (KP) equation, which are prototypes of what is called a soliton 
equation. For finite dimensional systems, the notion of integrability is clearer, i.e. it 
implies the existence of enough number of conserved quantities so that the initial value 
problem can be solved. The classic examples as Euler, Lagrange and Kovalevskaya 
tops belong to this category. We remark that the Toda equation is also an important 
dynamical system which is integrable either on finite or infinite lattices. 

And so what about the BBS in which we have just observed "solitons"? Are they 
really solitons that possess the above mentioned properties? Is BBS really integrable 
in some sense? Is it related to integrable systems known hitherto? Is there any good 
mathematical framework to analyse it? 

The aim of this review article is to give an introductory exposition on a variety of 
aspects of BBS elucidated in the last two decades, where all the above questions will be 
answered affirmatively. 

1.2. Overview of related mathematics 

It turns out that BBS originates in a quantum integrable system as well as in a classical 
integrable system. It is located at the very special point where the two systems meet 
by the procedures called crystallization and ultradiscretization, respectively. 

By quantum integrable systems, we mean those associated with the Yang-Baxter 
relation [6l [M]. Their symmetry is governed by quantum group Uq = Uq{g) meaning 
the g-deformation of the universal enveloping algebra U{q) of some affine Lie algebra 
g [TT| |33] . Typical examples are solvable lattice models in statistical mechanics such 
as the six- vertex model [6] for g = si2 and its generalizations. They are spin systems 
whose Boltzmann weights are continuous functions of q. The crystallization corresponds 
to taking the limit g — )■ 0, where the models are frozen to the ground state and its profile 
turns out to reproduce the BBS dynamics exactly. 

By classical integrable system, we mean here integrable difference equations such as 
discrete KP equation and (time-discretized versions of) Lotka-Volterra equation, Toda 
equation and so forth. These equations are already defined on lattices (discrete space- 
time), but their dynamical variables are yet continuous. The ultradiscretization is a 
procedure to transform these nonlinear evolution equations and their solutions into 
piecewise-linear forms. Leaving technical cautions aside, it is achieved by switching 
from the original variable a to A hj a = e~^^^ and taking the limit e +0. Being 
piecewise-linear, the resulting equations allow one to restrict the dynamical variables 
to a certain discrete set. In this way one reproduces evolution equations of BBS and 
obtains their solutions. 

Having the double, classical as well as quantum, origins of the integrability makes 
the study of BBS especially rich. One can import a variety of notions and techniques 
to understand and analyse BBS from the two theories. For instance from the theory of 
quantum integrable systems, we have ingredients like Yang-Baxter relation, quantum R 
matrices, commuting transfer matrices, Bethe ansatz, corner transfer matrices and so 
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forth [HI Eni [Ml SSI EH ES]- Similarly, the classical theory provides us with solitons, the 
inverse scattering method, tau functions, action-angle variables, isolevel set, spectral 
curves, linearization of flows etc [H IH El HSl [151 [23 [601 [ZS] • It turns out that they all 
survive the crystallization or the ultradiscretization rather miraculously. Moreover, they 
allow a systematic (Lie algebraic) generalization beyond the original BBS so that there 
are many kinds of balls or particles/anti-particles, boxes with capacity greater than one, 
and a family of commuting time evolutions etc. (Nonetheless, we will mainly focus on 
the basic type s[„+i case in this paper to be introductory.) 

Compared with traditional integrable systems, certainly a novel feature of BBS 
is that its dependent (or dynamical) variables have also been discretized. This fact 
indicates and actually has led to a fruitful connection to the realm of combinatorics. 
From a mathematical point of view, the crystallization and the ultradiscretization are 
both connected and actually have partly motivated the fascinating subjects known as 
crystal base [391 [10] iii the representation theory of quantum group, geometric crystals 
[7] as its geometric counterpart, and tropical geometry [591 [S3] in algebraic geometry. 
As the title of the article suggests, this review also contains elementary expositions and 
practical applications of these theories to BBS. 

Leaving the details to later sections, we present a rough schematic view of the 
relevant subjects as a summary. 

vertex model BBS integrable difference eq. 

quantum group crystal base theory geometric crystal 



1.3. Contents 

The main contents of each section follows: in §2.1 the basic notion of 

crystallization is illustrated with the simplest vertex model of ?7g(sl2). §2.2 is an 
exposition on sin+i crystal base theory, which is applied to describe the infinite BBS in 
§2.3. In §2.4, we briefly sketch various generalizations of BBS associated with affine Lie 
algebras. We remark that §2.2 and §2.3 are essential to study BBS. 

In §3.1 the physical background of Kerov-Kirillov-Reshetikhin (KKR) bijection in 
the Bethe ansatz is explained. The definition of KKR bijection for sl„+i crystal is given 
in §3.2 with the concrete algorithm. In §3.3 we state that KKR bijection linearizes the 
time evolution of BBS, which enables us to solve the initial value problem of BBS. 

In §4.1 the notions of tropicalization, ultradiscretization and min-plus algebra 
are introduced. In §4.2, two kinds of evolution equations for BBS are provided, 
corresponding to the "spatial" and the "soliton" descriptions of BBS. The equations 
are the ultradiscretization of known integrable discrete systems. The first description 
is studied in §4.4 and §5.2, and the second one is studied in §6.3. In §4.3, we briefly 
explain geometric crystal whose ultradiscretization gives the crystal structure. In §4.4 
the general solution for BBS is given by the ultradiscrete tau function in a similar way 



tropical geometry 




algebraic geometry 
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to many soliton equations. 

In §5.1 the basic features of periodic BBS is explained and its general solution is 
constructed via (modified) KKR bijection in §5.2. A remarkable feature is that the 
solution can be written in terms of tropical theta functions. In §5.3 we discuss more 
detail of periodic BBS from the viewpoint of torus decomposition of the isolevel set and 
fundamental periods of the time evolution. 

§6.1 is an introduction to tropical curve theory which is the latest mathematical 
object in this article. This theory is applied in §6.2 to solve the tropical periodic Toda 
lattice (trop-pToda). In §6.3 we show that the isolevel set of the periodic BBS is 
embedded in that of trop-pToda. This embedding bridges two different approaches to 
the periodic BBS by the trop-pToda and by the modified KKR bijection in §5.2, from 
the viewpoint of the tropical geometric description of the isolevel sets. 

We show a flow chart of the sections in this article: 
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where sections involving somewhat advanced or specialized topics are indicated by *. 

We did not intend to make the reference exhaustive. It is a moderate but sufficient 
supply for interested readers to proceed and find further references. 



2. BBS and crystals 

2.1. Crystallization: g — > of the vertex model 

For simplicity, we concentrate on the Uq{sl2) case in this subsection. Let us recall the 
six-vertex model and its fusion. Consider the two dimensional square lattice, where each 
edge is assigned with a local variable taking values in {1,2}. Around each vertex, we 
allow the following 6 configurations with the respective Boltzmann weights: 



1 2- 



2 2- 



■2 1- 



1 1- 



■2 2- 



12 12 12 

1 — q'^z q{l — z) q{l — z) z{l — q^) 1 — g^. 



(2.1) 



where z is called a spectral parameter. The other 10 configurations are assigned with 
Boltzmann weight. Let V = Cvi © Cv2. Then (12. ip is arranged in the quantum R 
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matrix R{z) G End(l^ ® V) as 

R{z) = a{z) Eii ® Eu + b{z) ^ Eu ® Ejj + c{z) (^zJ2 + J2) ® 



i<3 t>j 



(2.2) 



a{z) = 1 — q'^z, b{z) = q{l — z), c{z) = 1 — q^. 



Here the indices run over {1,2} and Eij is the matrix unit acting as EijVk = SjkVi. 
Schematically (12.21) is expressed as 



-^i )Eij (g)Eki, 



ijkl }^ ijkl I 

where the ^-dependence is exhibited. The Yang-Baxter equation 

R2^{z')R^^{z)R^2{z/z') = Ruiz/z')Rrsiz)R2,{z') 

12 3 

holds [6], where the indices signify the components in the tensor product as V 
on which the both sides act. It is depicted as 





(2.4) 

The R matrix R{z) is associated with the quantum affine algebra Ug = Uq{sl2) 
[TT| 133] . There is an algebra homomorphism A : Ug ^ Ug <^ Ug called coproduct, 
which enables one to construct the tensor product representation V from any two 
representations V and V'. Setting R{z) := PR{z) with P being the transposition of 
components, the quantum R matrix is characterized by the condition A{x)R = RA{x) 
for any x ^ Ug. The asymmetry between the last two in (12.11) is due to the special choice 
of the coproduct A that suits the limit g — that will be considered in what follows. 
(The precise form of A is not needed in this paper.) 

Starting from the six- vertex model, one can construct multi-state ("higher spin") 
solvable vertex models by the fusion procedure [17]. Let be the irreducible Ug 
module spanned by the m fold g-symmetric tensors in V®"^. It is a g-analogue of 
the spin y representation. Concretely, V\ = V and Vm with m > 2 is realized as 
the quotient l/®'"/^, where A = ^ ■ 1/®^' ® lmR{q-^) (g) V'®'"-2-i. it is easy to see 



Im R{q 



Keri?(g2 



(fi ® f 2 — qv2 ® vi). We take the base vector of Vm as 



^vf^'^ mod A, where Xj G Z>o and X1+X2 = m. The base will also be denoted by 



X2 



X 



(xi, X2) or by the sequence 1 ... 1 2 ... 2. Obviously dim Vrr. 



m 



1. The outcome 
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of the fusion procedure is the fusion R matrix R^"^'^\z) G End(K„ (g> Vi) given by 

J 



-y y^Vk, 



fc=l,2 



3 

x—r^y 
k 



^1 -g2x-i-)^ 



-.2X2 



(J,^) = (2,l), 
(J,A;) = (1,2), 



(2.5) 



(2.6) 



where y = (^1,^/2) is specified by the weight conservation (so called "ice rule") as 
yi = Xi + 6ij — 6ik. The RHS of (12. 6p is to be understood as unless this condition 
is satisfied. For m = 1, one has R^^'^\z) = R{z) and (12. 6 p reduces to (12.10 . For m = 2, 
(12. 6 p reads explicitly as follows: 

11^^11 12^^11 22^^12 12^^12 22^^22 



1 — q z 1 — q 



2 



q — q^z q^ — qz 



11- 



■12 12- 



■22 22- 



■22 11- 



■11 12- 



■12 



1^1 — q^)z (1 — q^)z 1 — q^z q^ — qz q — q^z . 



(2.7) 



Here we have suppressed z in the diagrams. (12. 7p is regarded as the (allowed) local 

configurations and their Boltzmann weights in a new vertex model where the horizontal 

and vertical edges take the 3 states {11,12,22} and the 2 states {1,2}, respectively. 

The weight conservation of R'^'^'^\z) means that the total number of letters 1 and 2 are 

preserved from NW to SE. 

Let us sketch how (12. 7p is obtained from (12. ip . The Yang-Baxter equation (12. 4p 

1 2 

with z' = zq^ shows that Imi?(g~^) C V ® V is preserved under the action of 

12 3 

Ri2,{zq^)R2z{z). Therefore its action on {V ® V) ® V can be restricted to V2 (8) Vi = 
[{V (g) V)/lmR{q-^)) (g) V. This yields the 2 x 1 fusion leading to R'^^''^\z). Similarly, 
can be deduced by restricting the composition {a{z) defined in (12. 2p ) 

i?i,„,+i(^g™-^)i?2,™+i(^g'"-=') ■ ■ ■ Rm.,m+iizq-'"^') 



a{zq"^-^)a{zq"'-^) ■ ■ ■ a{zq-"^+^) ^"^'^"^ 

to Vm'gVi. One can furthermore fuse R^^'^\z) along the other component of the tensor 
product in a completely parallel fashion. The result yields the quantum R matrix 
^(m,0(^2;) G End(Kn ® Vi). The R matrices so obtained again satisfy the Yang-Baxter 
equation in End(Vz ® ® Vk): 



R^-''\z')R'-^\z)R(r{z/z') = (./.')i?ir^(^)^: 



(l,m). 



23 



(2.9) 
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It is depicted as fl2.4p with the three hnes to be interpreted as representing Vi, Vm and 
Vfc. The quantum R matrix R'^'^'^\z) gives rise to a fusion vertex model on a square 
lattice in a similar manner to (12 .Sp and (12. 5p . The local variables on the horizontal and 
vertical edges are taken from Vm and V/, respectively. In terms of the linear operator 
:= PRi^'i){z) : Vm ^ Vi ^ Vi ^ Vm, which is also called an R matrix, the 
Yang-Baxter equation (12. 9p takes another familiar form: 

{R^'^'''\z') ® 1) (1 ® R^'^''\z)) {R^''"'\z/z') ® 1) 
= (1 ® i?^''") {z/z')) (^('''^^ (z) ® 1) (1 ® i?^'"''^) (z')) • 

It is (2;) rather than i?^''"'') (z) that will be directly related to the combinatorial or 
birational R introduced in the later sections. 

An important object in the vertex models is the (row to row) transfer matrix. For 
simplicity, we consider the basic case corresponding to R^'^'^'>(z). Then, the transfer 
matrix Tm{z) is defined by 

/ jl J2 jl \ 

Tm{z){v,,^---®v,J = J2Yl (^^'^4^^^'^4f Ux(i)]^^;i,®---®^fc„ (2.11) 

k, k2 kL 

where each ki runs over {1,2} and x^*-* does over the set of base of Vm labeled with 
{1...11, 1...12, . . . , 2. ..22}. The array of the vertex diagrams means the product of the 
corresponding Boltzmann weights (12.60 . We have assumed that the (horizontal) length 
of the lattice is L and employed the periodic boundary condition. The transfer matrix 
allows one to express the partition function of the model (for R^"^'^\z)) on N x L 
lattice with the periodic boundary condition as Z = TT{Tm{z)^). All the matrices 
Ti{z),T2{z), . . . act on the same space Vf^. Using (12. 9p . one can show that they form 
a commuting family: 

Tm{z)Ti{w)=Ti{w)Tm{z). (2.12) 

Now we are ready to discuss the main issue of the present section, namely, the 
crystallization limit g — j- 0. In (12. 7p . we see that only the 6 configurations in the left 3 
columns survive. In terms of the R matrix, the result may be stated that R^^'^^z) has 
the following action at g = 0: 

ll(g)lh^ 1®11, 12®lh^ 2® 11, 220 1^20 12, 

(2.13) 

11 ® 2 ^(1 ® 12), 120 2 ^ 2(1 ® 22), 22 ®2^2®22. 

Here {11, 12,22} and {1,2} are to be understood as labels of the bases of V2 and Vi, 
respectively. Apart from the factor z, (I2.13P provides a bijection between the two sets. 
The same feature can be checked easily for the general m case (12. 6p . (It is immediately 
seen for m = 1 by (12. ip .) The configurations having non- vanishing matrix elements 
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(Boltzmann weights) at g = are the following: 



2 



-1...12...2 1...12...2 ^— 1...12...2 

(0 < X2 < m) ^ (0 < X2 < m) . (2.14) 



The limiting Boltzmann weights are all 1 except the bottom right type, in which case 
it is z. The configurations fl2.14p determine a bijection between the data on the NW to 
SE, generalizing f l2.13p . In a physical terminology, the limit g — )■ corresponds to the 
low temperature limit, where crystallization takes place. Namely, spins are not allowed 
to thermally fiuctuate and are frozen to the ground state configuration determined from 
their choice at the boundary of the lattice. Here is an example of such a configuration: 

1 2 2 2 1 2 1 1 1 
111+111 + 112 + 122 + 222+122+222+122+112+111 

11112 12 2 2 
222+122+112 + 111 + 111 + 112+111+112+122+222 

2 2 2 112 111. (2.15) 

Regard such a configuration on 2 dimensional lattice as successive downward transfer of 
the horizontal array of spins on vertical edges. Then each step is a deterministic map 
corresponding to the crystallization of the transfer matrix (12. lip . The example fl2.15p 
corresponds to T-^il). This is an origin of the BBS time evolutions Ti, T2, . . .. The spins 
on the horizontal edges are "hidden variables" playing the role of carrier [73] . 

In the argument so far, one starts with g-dependent objects, e.g. fusion i?-matrices 
and transfer matrices, and then consider their crystallization g — )■ 0. In the subsequent 
sections, we explain how such procedures can be simplified and even more systematized 
by invoking the crystal theory of the quantum group Ug [391 SOI [28] . It provides a general 
framework to set up everything at g = from the outset. The labeling set of the bases of 
Vm and the quantum R matrix at g = will be formulated as crystal and combinatorial 
R, respectively. The power of z in fl2.13p is called the energy fl2.33p . which will also 
play an important role. The BBS and its generalizations will be constructed as the 
canonical dynamical systems associated with the crystalline vertex models. We remark 
that a similar approach to the BBS by the crystallization of the quantum Lotka-Volterra 
lattice has been undertaken in [251. 



2.2. Elements of crystal base theory 

The theory of crystal bases was founded by Kashiwara [39] BO] as a representation 
theory of quantum group Ug at q = 0. The notion of crystal is abstracted from the 
theory of crystal base [37] . In this subsection we give a brief description of crystals 
and combinatorial R which are basic ingredients in BBS. We remark that the notion 
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of crystal lattice in the theory of crystal bases (see [2S] for example) is omitted in this 
review, and our attention is focused on the notion of crystal. 

2.2.1. Definition of crystals. Let / be an index set. A crystal 5 is a set equipped with 
maps Ci, fi : B ^ B U {0} for i & I, satisfying certain axiom. In this article we are 
exclusively concerned with the semiregular case of [28| Def. 4.5.1]. Then the relevant 
axiom is the following: 

• For any b E B and i E I, there is > such that e"6 = 6 = 0, 

• eiO = fS = 0, 

• For bi,b2 G B, fibi = 62 is equivalent to 6^62 = ^i- 

Here we have omitted the items in the axiom involving the weight wt not used in this 
article. The and fi are called Kashiwara operators. They serve as g = analogues of 
Chevalley generators. For b E B, we set 

e,{b) = max{m\{e,r{b) ^ 0}, <^,(6) = max{m| (/,)'"(&) + 0}. (2.16) 

For our construction of BBS, we use the crystal Bi associated with sl^+i, where we 
take J = {0, 1, . . . , n} with n G Z>i. As a set Bi is given by 

Bi = {{xu . . .,Xn+i) G Z"+i I X, > 0, = " (2-17) 

The elements of Bi are also represented by Young tableaux. For each x = 
{xi, . . . ,Xn+i) G Bi we associate a one-row semistandard tableau of length / in which 
letter i appears Xi times. For instance let n = 2 and / = 2. Then the crystal 

B2 = {(2, 0, 0), (0, 2, 0), (0, 0, 2), (1, 1, 0), (1, 0, 1), (0, 1, 1)} (2.18) 

is also written as 



^2 = {1 1,22,33,12,13,23 }. (2.19) 



In what follows all indices of Xi,yi, . . . are interpreted in Z„,+i = Z/(n + 1)Z, namely 
Xi+n+i = Xi. The Bi is the labeling set of the bases of the Z-fold symmetric tensor 
representation (an example of so-called Kirillov-Reshetikhin modules) of f/g(s[„+i). 

For X = {xi, . . . , Xn+i) G Bi, let Cj, fi : Bi ^ BiU {0} (0 < z < n) be maps defined 

by 

ei{x) = {...,Xi + l, Xi+i - 1, . . .), /i(x) = (..., Xi - 1, Xi+i + 1, . . .), (2.20) 

if their images fall into Bi, or they are interpreted as otherwise. According to f l2.16p . 
the maps ei,ipi : Bi ^ Z [0 < i < n) are given by 



ei{x)=Xi+i, ipi{x)=Xi. 



(2.21) 
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For any crystals B, B' one can define their tensor product B ® B'. As a set it 
is a direct product B x B', but it also has a crystal structure. Any (x, y) E B x B' 
determines an element x ® y € -B ® B', and we understand x®0 = 0(8>y = 0. For 
X ® y E B ® B' the maps e^, y^j, gj, /j are given by 

ei{x ®y)= + {ei{y) - (2.22) 
(^j(x ®y)= i^iiy) + {<^i{x) - ei{y))+, (2.23) 

I CjX ® y if o9j(x) > £«(?/), 
e,(x®y) = <^ / ./ " (2.24) 

[X (g) Ciy if (y9j(x) < £^(2/), 

[xO/i?/ if V2i(x) < £»(?/), 

where (x)+ = max(x, 0). The tensor product defined in this way satisfies the axioms 
of the crystals. By repeated use of this construction one can define tensor products of 
more than two crystals, where the (co) associativity {B (g) B') ® B" = B ® {B' B") 
holds. In particular, this allows one to define the s[„+i crystal Bi-^ ® ■ ■ ■ ® Bi^ for any 
set of positive integers /i, . . . , /„. 

The crystals are represented by colored oriented graphs, known as crystal graphs. 
Let us show an example. 



B,: 





2 


B2: 


1 


1 




1 


2 




2 


2 


1 










1 






1 







(2.26) 



Here the arrows with index i represent the actions of /j. A tensor product of crystals is 
represented as follows. 



Let us show two more examples. 



2 ® 2 







1 




1 










2 




1 




2 




1 


2 



(2.27) 



Bi®B2: I 1 I 1 



° r^^ iii2ii — " r^® i2i2 



2|0 |2|2 
/ 



f2.28) 
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B2®Bi: I 1 I 1 I ® \l] 



2 2 ® 2 





1 






1 










1 




2 




1 


2 




2 



(2.29) 



2.2.2. Combinatorial R and its explicit formula. In general, two crystals B ® B' and 
B' ®B share a common crystal structure. The combinatorial R is the bijection between 
B ® B' and B' ® B that commutes with the actions of Kashiwara operators. It is 
a g — > limit of the quantum R matrix R{z) in §2.1. In other words it is a map 
R = Rbb' '■ B ^ B' B' B which satisfies the following relations: 



R{ei{x ® y)) = ei{R{x ® y)), R{fi{x ® y)) = fi{R{x ® y)). 



(2.30) 



In all the cases we consider in this article, the combinatorial R is uniquely determined 
by demanding the above conditions. 

By the definition, the inversion relation Rbb' ° Rb'b = Mb'^b holds. The simplest 
case is B = B', where the combinatorial R reduces to the identity map. As a non-trivial 
example, we find that the combinatorial R : B2®Bi Bi®B2 for SI2 is given by (12.131) 
by comparing the crystal graphs (I2.28P and (12.291) . modulo the power of z (which will 
be related to the energy function in (I2.33P ). For m general, R : Bm ^ Bi ^ Bi ^ Bm 
for is given by (12.141) in the notation of (I2.37p . More generally, we shall present a 
simple algorithm for R : Bi ® Bii ^ Bii ® Bi in §2.2.31 

For the sl^+i crystals there is a piecewise-linear formula for the combinatorial 
R. Given x = {xi, . . . ,Xn+i),y = {yi, . . . ,yn+i) e let x = {xi, . . . ,Xn+i), 

y={yi,..., yn+i) e ^""^^ be defined by 



Xi = Xi- Pi{x, y) + Pi_i(x, y), yi = yi + Pi{x, y) - Pi^i{x, y), 

(n+l k 
j=k j=l 



(2.31) 



Proposition 2.1 Given x E Bi,y G Bii define x, y E Z""*"-*^ by \2.31\) . Then: 
(i) All their elements are non-negative, hence x E Bi,y E Bf . 

(a) Define R : Bi ® Bii ^ Bi/ Bi by R{x ® y) = y ^ x. Then it is the combinatorial 
R for the sl„+i crystals, i.e. it satisfies the relations Ii2.30\) . 

One can prove it by showing the equivalence of the piecewise-linear formulas (I2.3ip with 
an algorithm for the R in §2.2.31 ([21 1 Prop. 4.1]). Another proof will be given, following 
an idea in ^91 Th. 4.28], as a consequence of the corresponding assertion in geometric 
crystals (Proposition I4.17p . 
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The formula of the combinatorial R fl2.31l) is characterized by the following relations 

Xi + Vi = yi + Xi, max(-Xj, -t/i+i) = max(-yi, -Xj+i), (2.32) 

with an extra constraint Y17=iixi ~ Xi) = Yl^=iiyi ~ Vi) ~ 0- '^^^ relations fl2.32p are 
consequence of ei{y<^x) = ei{x®y) and ipi{y^x) = ipi{x®y) which follow from ( I2.30p . 

A notion related to the combinatorial R is the energy function H : B ® B' ^ 
For the s[„+i crystals, : i?/ (g) i?// — Z is explicitly given by 

H{x®y) = Po{x,y) -mdjc{l,l'). (2.33) 

Even ii B = B' where we have R = Id, the energy function is not trivial and plays an 
important role in the theory of crystals and its applications [371 EH]- For the BBS, H 
will be used at Proposition 12.131 



2.2.3. Algorithm for combinatorial R. There is a simple way to calculate the image of 
the combinatorial R and the energy function without drawing the whole crystal graph, 
due to Nakayashiki and Yamada j6l] . We explain the algorithm along an example: 



1 


3 


3 


4 


7 




1 


3 


5 




1 


4 


7 


® 


1 


3 


3 


3 


5 



(2.34) 



Given the left hand side we can obtain the right hand side by using the following diagram: 




••• 



We suppose / > /' but to guess the algorithm in the case / < /' is easy. Represent 
X = (xi, . . . , Xn+i) G -B; by a pile of n + 1 boxes in which there are dots in the ith 
highest box. Do the similar for y = {yi, . . . ,yn+i) G -Bf and then juxtapose these piles 
of boxes. Repeat the following procedure (l)-(3) to obtain /' pairs of connected dots. 
(All the dots are unconnected initially.) (1) Choose any unconnected dot A in the right 
pile. (2) Look for its partner B in the left pile which is an unconnected dot in the 
lowest position but higher than that of A. If there is no such dots, B is chosen among 
unconnected dots in the lowest position. (We call the former case unwinding and the 
latter winding.) (3) Connect A and B. At the end we transfer all the unconnected dots 
from the left pile to the right one horizontally, yielding the piles for R{x ® y). The 
energy function f l2.33p is given by 

H{x ® y) = 7^(winding pairs). (2.35) 

For the above example fl2.34p we have H{x ®y) = Pq{x, y) — max(/, /') = 6 — 5 = 1. 

The algorithm for R and H will serve as the most substantial tool to check the 
examples in g2231 and 
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2.2.4- Yang-Baxter equation. The most important property of the combinatorial R is: 



Proposition 2.2 The following relation holds on B ® B' ® B" : 

{R ® ® R){R ® 1) = (1 ® ® 1)(1 (g) R). 



(2.36) 



The relation f l2.36p is known as the Yang-Baxter equation. We depict the relation 

R{x ® = y <g X hj 

y 

X- 



-X 



y 



(2.37) 



Example 2.3 By using the algorithm in ^2.2.3[ one can observe that the maps in the 
both sides of l\2. 36\) send an element in Bq^B^^^Bi to the same element in Bi®B^®Bq. 



2 


2 


3 


4 


5 


5 


® 


3 


3 


4 $ 


?) 6 



(i?® 1)(1 ® i?)( 2 2 3 



3 3 4 4 5 5 



6 ) 



2 2 5 





5 


® 


3 


3 


4 


4 


5 


6 ) 










3 


3 


4 


4 


5 


6 





l®R){R®l){l®R){ 223455033406 



2 2 3 


4 


5 


5 


® 


4 


® 3 3 


6 



'1®R){ 3 



2244550336 



3 


3 


4 


4 


5 


6 



It is also depicted as the following diagrams, where the lines represent crystals and their 



crossings stand for the combinatorial Rs. (See ^2.37\ ).) 
334 



223455 



334455 



223 



334456 



223455 
3 



225 



334 



336 



224455 



225 



334456 



(2.38) 



2. 3. Basic features of BBS 

In this section we introduce a one-dimensional cellular automaton associated with s[„+i 
crystals. For a more extensive presentation, see |7Uj . 
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2.3.1. States and time evolutions. For any positive integer L, we define a dynamical 
system on which generahzes the BBS in §1. In the system one may regard 

[~a~| G -Bi (a > 1) as a box of capacity one containing a ball with color a inside it, and 

p~| G Bi as an empty box of capacity one. We call our dynamical system an sin+i 
BBS. It is a cellular automaton equipped with a family of commuting time evolutions 
Ti, T2, . . . defined in the sequel. 

Bi ^ Bi — >■ Bi ® Bi be the combinatorial R and define Ri = 



Let R 

i-l 



Bi ^ 

L-i 



Id- 



Id Oi? O Id (g) ■ 



® Id {1 < i < L), which is a map from (_Bi)'^*~-^ 
I) Bi(g) {Bi)®^-\ Then Tl = Rl o ■ ■ ■ o R^ is a. map from 
^ to {Bi)®^ ® Bi. Given an arbitrary v ® hi ® ■ ■ ■ ® hi e Bi ® (Bi)'^^ let 



Bi ® to (5i 

Bi ® (5i 



7^(^; ® 61 ® 



hi 



® h'^ ® f '. It is depicted by 

hL^i 



-Vi 



b'l 



-V2- 



-Vl^2- 



-Vl^i- 



or simply by 



hi 



b'L-i 
hL 



(2.39) 



h[ h'^ 



v'. 



h'r 



We assume that the conditions 

L > 1, 

are satisfied in fl2.39p . and take 



6- =|X| for alH > 1, 



ui :-- 



Then we have v' = ui and the set {6'^, . . . , 6^} coincides with . . . 
the order of its elements gets shuffled. Under this setting let Ti : [Bi] 
El : (i?i)®^ — )■ Z be the maps given by 



Ti{hi ( 
Ei{hi 



hL) = h'l 



h' 



Li 



hL) = Y,{l-H{v,.i®h,)), 



(2.40) 



(2.41) 



(2.42) 

6l} as a set but 
^ (Bi)®^ and 



(2.43) 
(2.44) 



i=l 



where Vq = v. Call T; the /-th izme evolution and i?; the /-th energy. We note that every 
summand of f l2.44p vanishes for z ^ 1 because of H{ui ® [T] ) = 1, which ensures the 
convergence of the energy in the limit L — )■ 00. 

In what follows we often use the symbol ^ to indicate that its two sides are 
transformed to each other by the isomorphism (composition of combinatorial Rs) of 
crystals. For instance, the relation (12.371) is expressed as x (8> y ^ y (>>> x. This notation 
allows one to write the relation fl2.43p as a crystal version of "Lax equation" 



Mi ~ Ti{p) (g) ui, 



(2.45) 
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where p = 61 ® ■ ■ ■ ® ^l- An element in Bi is regarded as a carrier which can carry at 
most / balls. (The notion of carrier was introduced in [73] in the case of s[2-) In the 
carrier x = (xi, . . . , G Bi, the count of balls with label i{> 1) is Xj. The Ui fl2.42p 

corresponds to a vacant carrier. In f l2.39p . a carrier runs from left to right, changing 
itself as — t'l — )■ — > ■ ■ ■ • Although it is nothing but a repeated use of the algorithm 
in §2.2.31 with I' = 1, one can regard it as a successive loading/unloading processes of 
balls into/out of the carrier. 

To illustrate how the carrier works, as well as how the energy f l2.44p is evaluated, 
we show a few examples for f l2.39p . 

Example 2.4 Carriers with capacity 4, 3 and 2. Consider the state att = 3 in Example 
\2.11\ below. By the time evolution T4 it evolves into the state at t = A as follows. 



1 1 1 1 — |- 1 1 12— |— 1 122— |— 1 222— ^ 



2 1 
-2222— [■ 

1 



1 3 2 4 3 3 

■j— 1222— |— 1 1 22— |— 1 1 23— |- 1223— |— 1224— |— 1 234— |- 1 334 
2 2 2 1 3 2 2 



We added • to each vertex which scores +1 to the energy \2.44^ . Assuming that there 

are only 1 's to the right of the first row, there are no more scoring vertices. Hence we 

have £'4 = 9. At the vertex with •, no winding pairs occur in the algorithm in 

which makes the value of the energy function Ii2. 35\) zero. 

If the state were evolved by T3, then the diagram would be as follows. 



Ill— |— 112— 1-122— ^ 



2 2 1 

222 — |— 222- 
1 2 



-122 — |— 112 — 1-113 — |— 123 — |— 124 — 1-134 — |— 334 
2 2 1 3 2 1 



Hence E-^ = 8. In the same way, T2 would change the state as follows. 

3 2 4 3 



11^12^ 



2 2 2 

22 — |— 22 — |— 22 

1 2 2 



1 1 
12 

2 



— |— 11 — 1-23 — 1-24 — |- 

2 113 2 



34 



-33 



Hence E2 = 6. It is easy to see Ei = 3 and Ei = 9 for I > 4. 

Based on the Yang-Baxter equation f l2.36p . one can prove the commutativity of the 
time evolutions and the conservation of the energy. 

Proposition 2.5 [T71 Th.3.2] The following relations are satisfied 

TiTiip) = Ti,Ti{p), EiiTiip)) = Ei{p), (2.46) 

for any IJ'>1 and state p = hi® ■ ■ ■ ^hi. 

The /-th energy Ei is the conserved quantity associated with the time evolution T;. 
Under the conditions fl2.4ip we define Too just by taking a formal limit / — > 00 in fl2.45p . 
In fact, T;(p) = Too{p) holds if and only if / is greater than or equal to the maximum 
amplitude of the solitons contained in p. (This is due to Theorem 13.71 especially f l3.19p . 
See ( I2.47P for how to determine the amplitudes of solitons.) Remark that Ti serves as 
a shift operator which moves every ball to its right adjacent box. 
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Time evolutions of a state will be illustrated by drawing Tl{p),Ti~^^{p),Tl~^'^{p), . . . 



downwardly. For a state T/(p) 



we omit the symbol (8> and write it 



Hence a- denotes the value at site i and time t. In what follows we write 
a- = . (a dot) instead of a- = 1 for simplicity. Here is an example of the time evolution 
under Too: 

t=3 

t=4 

t=5 



.322554433.6 

322. . .5564433 

322. .5 654433. 



The special time evolution Too admits an elementary algorithm due to Takahashi 
generalizing the SI2 case in §1, which is a (non-local) description without a carrier. 

Proposition 2.6 Too = ■ ■ ■ Kn^i where Ka is an operator that works as follows, 

(i) Exchange the leftmost a with its nearest right [T], 

(a) Exchange the leftmost a among the rest of the [^'s with its nearest right p~| . 



(Hi) Repeat (ii) until all of the a 's are moved exactly once. 

K2K^K/^Kr,KQ to the t = 4 state in Example \2.15\ to 



Example 2.7 We apply 
obtain the t = 5 state. 

t=4 



.322. 
.322. 
.322. 
.322. 
. .223 



5564433. 



55.44336 

. .5443365 

. .5. .336544. . . 
. .5 654433. 



t=5 



.322. .5. 



,654433. 



A comparison of our formalism of BBS with that of the vertex models is summarized 
in the following table. 





vertex models 


BBS 


local states 


t/q-module 


crystal 


local interaction 


quantum R 


combinatorial R 


Ti 


transfer matrix 


time evolution 



2.3.2. Solitons. Now we define solitons in BBS. Intuitively, a pattern like ii...ii 
satisfying the condition ii > ■ ■ ■ > ii > 1 can be regarded as a soliton of amplitude 
/. It is denoted by [ii . . .ii]. In Examples 12.81 and 12.91 below, the sequence 554322 is a 
soliton of amplitude 6. 

Example 2.8 Time evolution by Ti with I > 6: 
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t=0 554322 

t=l 554322 

t=2 554322 

Example 2.9 Time evolution by T4; 

t=0 554322 

t=l 554322 

t=2 554322 



If well separated from the others, a soliton of amplitude / travels at a speed of min(/, k) 
under the time evolution T^. In particular we have the following for any 1-soliton state: 



Proposition 2.10 [T71 Lemma 4.1] If there exists only one soliton in the state, it travels 
at a speed o/min(/, k) under the time evolution where I is the amplitude of the soliton. 
And the value of the associated energy Ek is also given by min(/, k). 

Under the intuitive definition of solitons, one observes that the number of solitons of 
each amplitude may look changing during their scattering processes. We rather want to 
treat solitons as conserved quantities in BBS. For the purpose, we define mi by 

El = min(/, k)mk or equivalently mi = —Ei^i + 2Ei — -E^+i, (2.47) 
fc>i 

where £"0 = is understood. Since Eis are conserved quantities, so are m^s. Then one 
can interpret mi as the number of solitons of amplitude /. In view of Proposition I2.10[ 
it is consistent with the previous definition of solitons in the case where solitons are well 
separated. By the definition, any state p is an A^-soliton state, where is determined 
by iV = Ei{p). We note that according to (2.42), Ei{p) equals the number of adjacent 

pairs [T]®[^ with a > 1 appearing in p, whereas -E'oo(p) is the total number of balls. 



Example 2.11 A three body scattering process under T^o- From Example \2.4\ we find 



att = 3 that mi = —Ei^i + 2Ei — Ei^i = 1 for / = 2, 3, 4 and mi = in the other cases. 

t=0 2222 332.. 43 

t=l 2222 332.43 

t=2 2222... 33243 

t=3 2222.. 32433 

t=4 222.322433 

t=5 22.. 3224332 

t=6 22. . .322.4332 

t=7 22 322.. 4332 

t=8 22 322... 4332 



We will see that the nonlinear time evolutions of BBS are transformed into linear ones 
on the rigged configurations. See Example \3.S[ 
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We remark that, besides m;s, our BBS has additional conserved quantities fl3.20p . By 
using the crystals for anti-symmetric representations, one can show that the color degrees 
of freedom for any state of the BBS can be transformed into a "word" which does not 
change under any [69] . 

2.3.3. Scattering rules. The scattering of solitons in our BBS consists of the exchange 
of their internal degrees of freedom and the phase shifts. Although it is possible to treat 
general many-body scattering processes, we devote ourselves to the case of two-body 
scatterings for simplicity. In what follows we assume that the time evolution is given 
by Too. 

Example 2.12 A scattering process of two solitons with amplitudes I = 6 and I' = 3 in 
5(5 BBS. 

t=0 554322 422 

t=l 554322 422 

t=2 554322... 422 

t=3 5543.. 42222 

t=4 553 442222 

t=5 553 442222 

t=6 553 442222 

t=7 553 442222... 

Suppose at time t = 0, the state bears two solitons [ii ■ ■ ■ ii] and [ji' . . .ji]. Denote 
this two-soliton state by [ii . . . ii]x x [j;/ . . .ji]y, where x and y are the positions of their 
leftmost letters. For instance, we have [554322] 1 x [422] ig in Example 12.121 at t = 0. 
We assume / > /' and x <^ y. Then the former catches up with the latter and they 
eventually collide. Before the collision these solitons travel at speeds of / and /' cells 
per unit time respectively, so at time t we have [ii . . . ii]x+it x [jv ■ ■ ■ ji]y+i't- After the 
colhsion we have two-soliton state [j;/ . . ■ ji]y+i't~s x [ii . . .ii]x+it+5 where 5 is a phase 
shift. (This fact is also a part of the statements in the forthcoming Proposition 12.131 ) 
By the collision the larger soliton gets pushed forward and the smaller soliton pulled 
backward by an amount of 6 cells. 

The exchange of the internal degrees of freedom occurring here is governed by the 
combinatorial R, and the phase shift is essentially given by the energy function H. 

Proposition 2.13 [17, Th.4.6] Any collisions of two solitons asymptotically break up 
into two solitons. Let the two-soliton state be [ii . . .ii]x+it x [ji> . . . ji]j^+;'t at time t well 
before the collision, and [jii . . ■ ii\y+vt~5 x [ii . . .ii\x+it+5 for t after the collision. Then 
the phase shift 6 is given by 



6 = H 



h]® \Ji\-\Ji'\ ) +1'- (2.48) 
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And the exchange of the internal degrees of freedom is described by the combinatorial R 
for the sin crystal^ 



R ii ■■■ ii (8) ji ■■■ ji' 



Ji 



Ji' 



(2.49) 



The phase shift can be computed by f l2.33p or fl2.35p . It is always positive and take 
values between /' and 21'. 

Example 2.14 The scattering process in Example \2.12\ is expressed as [554322]i+et x 
[422]i6+3t —7- [553]ii+3t X [442222]e+6t- This is described by the combinatorial R for the 
crystals. By using the algorithm in ^2.2.31 it can be computed as 



•••• 



The value of the energy function (= # (winding pairs)) is 2 and we observe that the 
phase shift is given by 6 = 2 + 3 = 5. 

Example 2.15 Scattering processes of three solitons. There are three solitons [554322], 
[433], [6] at time t = 0, and again three solitons [3], [522], [654433] at time t = 8. 

t=0 554322 433 6 

t=l 554322... 433 6 

t=2 554322433 6 

t=3 322554433.6 

t=4 322. . .5564433 

t=5 322.. 5 654433 

t=6 3225 654433 

t=7 3522 654433 

t=8 3.. 522 654433 



t=0 
t=l 
t=2 
t=3 
t=4 
t=5 
t=6 
t=7 
t=8 



554322 . 



.554322 



433. .6 

4336 

554322 4633 

554322 4. .633 

5543242. . . 633 

3.554422633 

3 522654433 

3 522. . .654433 

3 522 654433 



:|: The elements of crystals which appeared in the formulas in Proposition 12.131 have no letter "1" . 
Hence the combinatorial R and the energy function used there are regarded as those for sin crystals by 
reducing the value of all letters in the Young tableaux by 1. 
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Both processes have the same kinds of soUtons at t = 0. The orders of collisions 
occurring in subsequent times are different. The fact that the outcomes at t = 8 share 
a common soliton content can be viewed as a consequence of the Yang-Baxter equation 
/ fO^) . See ExamplelKM 

Due to the commutativity of the time evolutions (Proposition 12. 5p . the scattering rule 
remains unchanged when T^o is replace by T^ for any k > V . In this case we have the two 
body scattering [ii . . . ii]x+m\n(i,k)t x [jv ■ ■ ■ ji\y+vt [jv ■ ■ ■ ji\y+vt~5 x [z/ . . . n]x+mm{«,fc)t+5- 
Here the phase shift and the exchange of internal degrees of freedom are still given by 
(ElHl) and ([239]). 



2.3.4- symmetry. So far we have not mentioned the role of the Kashiwara operators 
acting on the states of BBS. Their significance is recognized as the s[„ symmetry in the 
system. Let p be a state of BBS. Suppose Cip ^ for some i G {2, . . . ,n}. Then we 
have 

T, (cip) = em (p) , El {e,p) = Ei (p) , (2.50) 

for any / [T7]. The first relation is a manifestation of the s[„ symmetry of the time 
evolution. Due to (I2.47P and the second relation in (I2.50p . this transformation does not 
change the amplitudes of solitons but alters their internal labels. 

A conserved quantity associated with the s[„ symmetry is defined as follows. Given 



a state p = ai ® ■■■ ^ cll ^ let Wi ... Wk he the word which is obtained from oi . . . by 
ignoring every 1. Denote by P{p) the P-symbol P{wk ■ ■ -Wi) and Q{p) the Q-symbol 
Q{wk . . .wi) (semi-standard Young tableaux) obtained from the opposite word Wk . . .wi 
by the Robinson-Schensted-Knuth correspondence [18]. Explicitly the P-symbol is 
defined as 

P{wk . . .wi) = Wk (wk-i ^ (. . . (ws iu;2 wi)) . . .)) 

where — >■ implies the column insertion and ^ does the row insertion [18]. The Q-symbol 
is the standard tableau consisting of {1, . . . , A;} that records the growth history of the 
P-symbol. We have the following: 

Proposition 2.16 [161 Th.3.1] The P-symbol P{p) is a conserved quantity of BBS, i.e. 
P{Ti{p)) is independent oft for any I. 

We note that the time evolution of the BBS is attributed to the dynamics of Q- 
symbol [161 Th.5.1]. 

For instance consider Example 12.151 The opposite words wio . . .wi are 6334223455 
for t = and 3344562253 for t = 8. Both words share a common P-symbol 

223455 
P{TUp))= 334 . 
6 
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The growth pattern of the pair (P, Q) by the successive row insertions for the word 
6334223455 looks as follows: 

6 1 3 1 33 13 334 134 
6 2 ' 6 2 ' 6 2 ' 

234 134 224 134 223 134 
3 2 , 33 26 , 334 267 , 
6 5 6 5 6 5 

2234 1348 22345 13489 223455 13489X 
334 267 , 334 267 , 334 267 
6 5 6 5 6 5 

Here X denotes 10. 

In the case of sl2 it is known that one can introduce another P-symbol whose shape, 
which represents the hst of amplitudes of solitons, is a conserved quantity [21 [80] . 



2.4- Various generalizations 

The BBS has been generalized extensively. Here we present a few prototype examples. 

• Generalizations in sin+i case. The original BBS consists of boxes with capacity 1 only, 
which corresponds to the fact that the states belong to ■ ■ ■ (8) Pi ® Pi (X> ■ ■ ■ . A natural 
generalization is to replace it with ■ ■ ■ (g) P^. ® Pfc^+i ■ ■ ■ (8) ■ ■ ■ . The commuting family 
of time evolutions {Ti} are defined in the same way as before, where the vertical lines 
in (12.401) now represent Pfc-S. The resulting dynamical system is BBS involving a box 
with capacity ki at the site i [2ll [78]. The basic features of the system, e.g. solitons, 
scattering rules, conserved quantities, linearization scheme, etc. remain the same as the 
capacity one case. See [2T|, [501 [M]- The BBS with a periodic boundary condition will 
be treated in §5. 

One can further use crystals other than the family {Bi}. Examples of such kind 
are a BBS with reflecting end [53] and a BBS associated with anti-symmetric tensor 
representations of Ug{5ln+i) [82] . 

• grt-automaton. Similarly to the stn+i(= An^) case, integrable cellular automata 
associated with the non-exceptional affine Lie algebra Qn = Bn \ Cn \ Dn\ A^2n-ii 
A\J^ and P„_,_i have been constructed [2l] and the soliton scattering rule determined 
[23] . The dynamics allows a neat description in terms of particles and anti-particles 
that undergo pair creations and annihilations [22] . The BBS turns out to be the special 
case in which no anti-particle is present. Let us demonstrate the ^4^^ case. Each 
local state takes values in {1,2,3,4,1,2,3,4}, where 2,3,4 are anti-particles of 2,3,4, 
respectively. As in the BBS, 1 represents an empty box whereas 1 plays the role of 
particle & anti-particle bound state. The prototype time evolution Too is given by 

Too = K2K,K,KiK-sK-2, (2.51) 

where each Ka is defined by the following algorithm (We understand 2 = 2 etc.): 
(i) Replace each 1 by a pair a, a within a box. 
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(ii) Move the leftmost a (if any) to the nearest right box which is empty or containing 
just a. (Boxes involving the pair a, a are prohibited as the destination.) 

(iii) Repeat (ii) until all of a's are moved exactly once. 

(iv) Replace the pair a, a within a box (if any) by 1. 

When anti-particles are absent, (i) and (iv) become void and the algorithm reduces 
to the one for BBS in Proposition 12.61 

Example 2.17 D^^ -automaton. We write ., a, b, c, d for 1, 1, 2, 3, 4, respectively. 



t=0 ..b22 3.. 

t=l b22..3. 



t=2 b223 

t=3 3c32. 



t=4 3..C32 

t=5 3 c32 

t=6 3 c32. 



t=2 b223 

a23 

32a 

324d . . . 

32. a. . . 

23c3. . 

t=3 3c32. 



A soliton is a consecutive array of the form 2'^2 3'^s4n/i'^i^'^3 2''2 ^ where z/j 's are nonnegative 
integers such that ugu^ = 0. [Left]: Successive time evolutions under T^, where pair 
annihilation/ creation b2 — > c3 takes place in the scattering. [Right]: The intermediate 
states between t = 2 and t = 3 corresponding to Ii2. 51\) . where the procedures (i)-(iv) 
can be checked. 

In general, it is expected that so-called Kirillov-Reshetikhin module has the crystal 
base [62] and one can use its crystal to construct the corresponding generalization of 
BBS. 

• Supersymmetric case. The supersymmetric automaton given by the crystal for the 
super Lie algebra A{m,n) was introduced in [26]. We have fermionic balls labeled by 
m-|-2,...,m + n + l, besides the empty boxes labeled by 1 and the {bosonic) balls labeled 
by 2, . . . , m + 1. The time evolution rule is the same as that for the sl„+m+i-automaton 
in Proposition 12.61 except for the step to move a ball with a fermionic label a. For a 
fermionic label a, we replace (iii) of Proposition 12.61 with 

(iii') Exchange the leftmost a among the rest of the a's with its nearest right 1 if this a 
has not been overtaken by the previously moved a. 

This rule denotes that each soliton can contain at most one fermionic ball of each label. 



Example 2.18 A{1,1)- automaton. The constraint in (iii') with a = 3 is relevant in 
the steps from t = 1 to t = 2 and from t = 2 to t = 3. 

t=0 ...322.. 3 

t=l 3223 

t=2 3322 

t=3 3.. 322 

t=4 3 322 
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3. Bethe ansatz approach 
3.1. Introduction 

The Kerov-Kirillov-Reshetikhin (KKR) bijection [HI US] is a one to one correspondence 

{rigged configurations} {highest paths}. (3.1) 

It originates in Bethe's consideration on the completeness of the Bethe ansatz under the 
string hypothesis [8] . We shall explain (13. ip after a brief exposition on the background 
along the simplest example from sl2- 

Consider the spin | Heisenberg chain with the Hamiltonian acting on (C^)*^^: 

L 

-n = + ^M+i + ^k^Ui - !)• (3-2) 

k=l 

Here is a Pauli matrix acting on the A;-th site and the periodic boundary condition 
al^^i = a I is assumed. The model possesses the (global) 5(2 symmetry in the sense that 
(j° := Ylik satisfies the defining relations of si-z, and [cr", %] = 0. Let = C^i © Cv2, 
where vi and V2 are regarded as spin up and down local states, respectively. As a 
consequence of the SI2 symmetry, the Hamiltonian T-L preserves the number of down 
(hence up as well) spins, so one may concentrate on a subspace Wr with r down spins 
and L — r up spins. The diagonalization of "H is done by the Bethe ansatz [8] . It reduces 
the task to finding the solutions of the Bethe equation (r < L/2): 

= ^fl^l^ltl^ ij = l,...,r). (3.3) 

In term of the Bethe roots {ui, . . . , Ur}, one can construct the eigenvector \ui, . . . , Ur) G 
Wr, called Bethe vector, of "H whose eigenvalue is given by YTj=i known that 

the Heisenberg Hamiltonian is contained in the commuting transfer matrices {Tm{z)}m,>i 
(I21T|1 with g = 1 as Ti{z) = Ti(l)(l + const(^ -l)n + ---) (cf. P §10.14]). Thus the 
Bethe equation (13.31) is actually relevant to their joint spectrum and therefore to the 
"diagonalization" of the commuting time evolutions in the BBS (although the latter 
corresponds to g = rather than g = 1). 

Back to (13. 3p . the variety of eigenvalues is provided by the variety of solutions to 
the Bethe equation. Thus a basic question arises; how many solutions should there be 
for the completeness of the Bethe ansatz ? The answer is (^) — The decrement 

from dim Wr = (^) is due to the fact that the Bethe vectors are by construction highest 
weight vectors annihilated by the 3(2 raising operator. Namely, it has the property 
a~^\ui, . . . , Ur) = with = {a^ + i(T^)/2 by construction [HI [8]. By virtue of the SI2 

§ In this argument, independence of the associated Bethe vectors has not been taken into account, 
and all the Bethe roots are supposed to be finite. 
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symmetry, the other eigenvectors can be produced by applying the lowering operator 
cr~ = [a^ — icr^)/2 successively. Thus one should be content with capturing all the 
highest weight vectors as Bethe vectors. 

Let us observe an example L = 6, r = 3. There are certainly (g) ~ (2) ~ ^ solutions 
as given below. 



0.8585 



-2.0175i 



-0.8585 



2.0175i 



Re Uj 



Im Uj 



0.4718 
-1.0006i 



0.4718 
+1.00061 



—1 



-0.9436 



-0.4718 
-1.0006i 



0.9436 



-0.4718 
-1. 00061 









n 









1 



Here each Bethe root Uj is depicted as •. Within each solution, they are grouped 
into strings. A string is an array of w's which is symmetric with respect to the real 
axis and equidistant of difference 2i with possibly "negligible" distortions. Strings 
consisting of k w's are called /c-strings. In the top left (right) solution, there are three 
1-strings (one 3-string). The three solutions in the bottom line consist of a 1-string 
and a 2-string with different real parts (called centers). These features are conveniently 
symbolized in a Young diagram (called configuration) where each row is attached with 
a nonnegative integer (called rigging) as shown in the figure. They are examples of 
rigged configurations. Each row including the rigging signifies the length and the center 
of the string encoded as an integer. They are to obey a certain selection rule that will 
be specified later in a more general setting. (See fl3.12p . The way to find the rigging 
will also be explained in section 13.21 ) To summarize so far, rigged configurations are 
combinatorial analogue of the pattern of Bethe roots under the string hypothesis. 

Let us turn to the RHS of (13.11) . Bethe vectors have the form \ui,...,Ur) = 
Yl (S> ■ ■ ■ ® f G Wr, where the sum runs over ii, . . . ,iL G {1,2} such that 

. . . , ii} = L—r and #2(^1, • • • , ii} = Highest paths are combinatorial analogue 
of them represented as the sequence zi, . . . , G {1, 2}^ satisfying the same condition 
as above and 

#i{?i, . . . , ^ J > #2{?i, . . . , for 1 < A; < L. (3.4) 

This is a remnant of the highest condition a^\ui, . . . = 0. There are (^) — [j.^^) 
highest paths as expected. In our example L = 6,r = 3, the highest paths and the 
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corresponding rigged configurations in fl3.ip are given as follows: 





< — ^121212 




]0 < — ^111222 







10 



i — ^121122 



Ul 



10 



i — ^112122 



Q2 



10 



i — >112212. 



(3.5) 



This is an example of the KKR bijection. The arrows — )■ here, or equivalently the map 
in f l3.ip . is a combinatorial analogue of the Bethe ansatz which produces Bethe 
vectors from Bethe roots as {ui, . . . , Uj.} ^ \ui, . . . , Ur). 

The (vague) claim that any solution of the Bethe equation can be described as a 
collection of strings is called string hypothesis. It is known that string hypothesis is not 
always valid literally (see for example [13]). Nevertheless, as we shall illustrate below 
for sln+i case, one can define rigged configurations and highest paths and establish their 
bijective correspondence mathematically. 



3.2. KKR bijection 



Let Bi be the sln+i crystal (I2.17P with / = 1. For simplicity we shall exclusively 
consider the crystal of the form Bf^ (L G Z>i) and call its elements as paths. For a 
Young diagram A with |A| = L and depth at most n + l, elements of the set 



V+{L, X) = {pe Bf^ \eip = (1 < i < n), wtp = A} 



(3.6) 



are called highest paths with weight A. Setting p 
(^iP = (1 < i < n) is concretely described as 



, the highest condition 
(3.7) 



...,ik}> #2{«i, • • • , 4} > ■ ■ ■ > #n+i{^i, • • • , 4} for 1 < A; < L, 

which is a generalization of (13. 4p . By the condition wtp = A we mean 

#4zi,...,Zi} = A, {l<a<n + l). (3.8) 

Let us proceed to the definition of the rigged configurations. Let . . . , /x*^"^ 

be an (r2 + l)-tuple of Young diagrams. We will always take /i'-'^^ = (1'^) {L E Z>i) which 
will match the choice of the crystal Bf^ in (13.60 . Denote by mj"'* the number of length 
j rows in /i^"-* and introduce the following: 



(a) 
P)' 



(a) 



q- ' --zqy +q- ■ ' {I < a < n), 



(3.9) 
(3.10) 



fe>i 



By the definition m^^ = Ldj^i and q^^' = L for j > 1. In general qj""' is the number of 



-,(0) 



(a) 



cells in the left j columns of fx^^\ The integer pj is called a vacancy and will play an 
important role in what follows. 
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An (n + l)-tuple of Young diagrams (/i'-"-', . . . , '-"■'') is a configuration if p^"^ > for 



any 1 < a < n and j G Z>i such that m^'^^ > 
rectangle constituting /i*^"^^) will be referred to as a block. 

Given a configuration . . . we attach a rigging Jj'^]^ G Z>o to every row 

in /x^*^-* except = (1^) as follows (shown for a block (a, j)). 

(^(a)^ J{a)) 



Such a pair (a, j) (i.e. the mj"'' x j 



m 



(a) 



J 



(a) 



t{«) 



(a) 



We group the rigging as J = ( J(^\ . . . , J^")) where J(") = (^^ )i>i,i<a<m 



(3.11) 

(a) is 



the one attached to A configuration ((1^), . . . j/U*-"-*) attached with a rigging 
J = (J(^), . . . , J(")) will be denoted by (/i, J)l with /i = (/i^^), 
is a rigged configuration if the condition 



< jjf < Jf,^ < ■ ■ ■ < /"^ 

— J.J^ — J.^ — — j^m 



M < Pj 



(a) 



We say (/x, J)^ 
(3.12) 



is satisfied for all the blocks (a, j)0. 



Example 3.1 We list all the rigged configurations having the configuration ({1^),^^^^) 
with n^^^ = (2, 1, 1). To save the space, only J^^^) part is given. 



□ 

2 




□ 

1 2 




□ 

2 2 




□ 

1 2 
1 



□ 

2 2 
1 



(3.13) 



For later convenience, we have exhibited the vacancy p^2^ = and p^^^ = 2 to the left of 
the relevant blocks. 

Example 3.2 An n = 3 example. Again, vacancies, e.g. p'^^ = 5, are exhibited. 



(0) 



Q14 



(1) 



(3) 









□ 1 












2 











3 





1 ono 



(3.14) 



II This condition is known actually to ensure that pf^ > for all j G Z>i. 
^ Equivalently, one may only impose < Jj"'-} , ■ ■ ■ , J^"^ ,„) < p^""^ and identify their permutations. 
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A weight of a rigged configuration J) ^ is the Young diagram A = (Ai, . . . , A„+i) 
specified by 

l/i^'^)] = Aa+i + Aa+2 + --- + A„+i (0<a<n). (3.15) 

We write it as wt((/i, J)l) = X, which is actually dependent only on the configuration. 
The inequality Ai > • ■ ■ > A„+i > is guaranteed by the condition p^{= — 
2|/i(")| + l/i^'^+^'l) > for 1 < a < n. Note also that |A| = L. Let RC(L, A) be the set of 
rigged configurations of weight A. 

Theorem 3.3 For any L G Z>i and a Young diagram A with |A| = L, there is a 
bijection 

RC(L,A) ^ V+{L,X). (3.16) 

</- 

The original KKR bijection [HI |13] is the one between rigged configurations and 
Littlewood- Richardson tableaux. Its ultimate generalization for type s[„+i is available 
in |l5l [65]. In the simple setting of this paper, the Littlewood- Richardson tableaux 
are in one-to-one correspondence with highest paths via the Robinson-Schensted 
correspondence [T8] . 

We regard a rigged configuration (/i, J)l as a multiset of strings. A string 
corresponds to a row in (13. lip . It is a triple (a,j, jj^^^) consisting of color a, length 
j and rigging jj'^^ A string is singular if jj'^ = pf'\ namely if the rigging attains 
the allowed maximum in (I3.12p . We regard the highest path p 



_ as a 

word iii2 ■ ■ - II G {1, . . . + 1}^. (The Littlewood-Richardson tableau mentioned in 
the above is the Q-symbol [18] of this word.) 

For simplicity, we first explain the algorithm for for n = 1 case. Even in 
this case, it may look formidably complicated at first glance. However, it is a very well- 
designed algorithm, and the readers will be impressed and get familiarized with it pretty 
well by working out a few examples. The nS^ and the vacancy p^^ will be denoted by 
nij and pj. Thus the definition (13. 9p becomes pj = L — 2 J2k ^^^U^ k)mk. It is useful to 
remember it as pj = L — 2(number of cells in the left j columns in the Young diagram). 

Algorithm of (p for n = 1. 

Given a highest path ii...iL G {1,2}'^, we construct the rigged configuration 
(j){ii . . Al) = (/i, J)l inductively with respect to L. When L = 0, we understand that 
0(-) is an empty Young diagram. Suppose that 0(ii . . .i^) = (/i, J)l has been obtained. 
We are to construct (/x', J')l+i = <P{ii ■ ■ ■ ^l^l+i) from (//, J)l and iL+i G {1,2}. 

Case i^+i = 1. One has (/i', J')l+i = (z^, J)l+i, which means that no change should 
be made in the length and rigging of the strings. (By the definition, their vacancies pj 
increase uniformly by one.) 

Case ii+i = 2. (a) If there is no singular string in (/i, J)l, just additionally create 
a length 1 singular string with respect to the new configuration. (Its rigging is therefore 
L + 1 — 2 min(l, A;)(mfc + 6ki)-) (b) If there exist singular strings, pick a longest 
one among them and let i be its length. (Any choice is OK when there are more than 
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one such strings.) Then (/i', J')l+i is obtained by extending the string to length i + 1 
and making it singular with respect to the new configuration. (Its rigging is therefore 
L + 1 — 2 mm{i + 1, k){mk — Sk/ + 6k/+i).) In either case of (a) and (b), keep the 
other strings unchanged. 

Algorithm of <p~^ for n = 1. 

Given a rigged configuration {fL,J)L, we construct a highest path = 
(f)~^{{fi, J)l) inductively with respect to L. We are to determine ii G {1,2} and 
in', J')l-i such that J)l) = J')l-i) ii- 

If (/i, J)l contains no singular string, then 2l = 1 and (/i', J')l~i = (/W, J)l~i- The 
latter means no change should be made in any string. (By the definition, their vacancies 
Pj decrease uniformly by one.) If {fi,J)L contains singular strings, then ii = 2. Pick 
a shortest singular string and let i be its length. (Any choice is OK when there are 
more than one such strings.) Then (/i', J')l-i is obtained by shortening the string to 
length i — 1 and making it singular with respect to the new configuration. (Its rigging 
is therefore L — 1 — 2^^min(£ — 1, k){mk + 5k,i-i — 5k/).) The other strings are kept 
unchanged. 



Example 3.4 For the rigged configurations in Example \3.1\ the algorithm of cf)^^ 
proceeds along the arrows. The algorithm of (p proceeds backward. To save the space, L 
is given in the first line. 













7 




6 






5 




4 




3 




2 




1 


lm(j) ^ 



2 




no 







1 




1 

















ino 


h 


ODO 


h 


3 1^ 


12121122, 



2 




no 
1 






1 




1 
1 








ll 


3J 




200 




ino 


h 


ODO 




3 1^ 


12112122, 



2 




no 

2 





imo 

300 




1» 


ll 


3J 




200 


h 


ino 


h 


ODO 


h 


3 1^ 


12112212, 



2 




no 

1 

1 




1 




1 
1 
1 






ll 


31 




2ni 


h 


ini 


h 







3 1^ 


11212122, 



2 




no 

2 
1 




imo 
stfi 




1» 




3} 




201 


h 


ini 







h 


3 1^ 


11212212, 



2 




no 

2 
2 


4 


I 






1» 


imo 




omo 


h 


ini 


h 





h 


3 1^ 


11221212. 



Aoie that one should keep updating the vacancies with L. 

Now we proceed to the n general case. The basic idea is to apply the 
removal/ addition procedure for the n = 1 case recursively in the direction of color. 
Algorithm of for general n. 

Given a highest path ii . . .i^, we construct the rigged configuration (f){ii . . .Il) = 
{fi,J)L inductively with respect to L. When L = 0, we understand that </>(■) is the 
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array of empty Young diagrams. Suppose that (f>{ii . . .il) = {h,J)l has been obtained. 
Denote z^+i G {1, . . . , n + 1} simply by d. We are to construct (/i', J') l+i = </)(zi . . . iid) 
from (yU, J)l and c?. \i d = 1, then (/z', = (/i, which means that no change 

should be made on any string. (By the definition (13. 9p . the vacancies p^^^ increase by 
5ai-) Suppose d> 2. 

(i) Set i^'^^ = oo. For c = d — 1, d — 2, . . . ,1 in this order, proceed as follows. Find the 
color c singular string whose length i^"^^ is largest within the condition i^'^^ < 
If there are more than one such strings, pick any one of them. If there is no such 
string with color c, set i^'^^ = 0. Denote these selected strings by [c,£^'^\ ji^^^ with 
c = d — l,d — 2, . . . ,1, where it is actually void when i^'^'' = 0. 

(ii) Replace the selected string [c,i^'^\ by (c, f*^'^-' + 1, ji*^^) for all c = d — l,d — 

(c) 

2, . . . , 1 leaving the other strings unchanged. Here the new rigging J; is to be 
chosen so that the extended string (c, f*^'^-' + 1, ji'^^) becomes singular with respect 
to the resulting new rigged configuration (yu', J')l+i. 

The algorithm is known to be well-defined and the resulting object gives the sought 
rigged configuration (/i', J')l+i = 0(^i • • -iLd). 
Algorithm of <p~^ for general n. 

Given a rigged configuration {fi,J)L, we construct a highest path ii...iL = 
0~^((/i, J)l) inductively with respect to L. We are to determine d{= ii) G {1, . . . , n + 1} 
and (/i', J')l~i such that 0~^((/i, J)l) = (t)~^{{fi', J')l~i) d. 

ii) Set = 1. For c = 1,2,... , n in this order, proceed as follows until stopped. 
Find the color c singular string whose length i^^"^ is smallest within the condition 
< £(c) jf are more than one such strings, pick any one of them. If there 
is no such string with color c, set d = c and stop. If c = n and such a color n string 
still exists, set d = n + 1 and stop. Denote these selected strings by {c,^'^^\ j'f^^ 
with c = 1, 2, . . . , c? — 1. 

{ii) Replace the selected string (c, j'f^) by (c, l^''^ - 1, J^"^) for all c = 1, 2, . . . , - 1 
leaving the other strings unchanged. When f'^'^^ = 1, this means that the length 
one string is to be eliminated. The new rigging ji^^ is to be chosen so that the 
shortened string (c, f*^'^-' — 1, ji'^'') becomes singular in the new data (/i', J')l-i- 

For an empty rigged configuration, we understand that 0~^((0, 0)^) = 

L 

0~^((0, 0)i_i) 1 = ■■■ = 11 . . . 1. The algorithm is known to be well-defined and 
ends up with the empty rigged configuration at L = 0. The resulting sequence gives the 
sought highest path ii . . Al = (t)^^{{fi, J)l)- 

Example 3.5 The algorithm of (p^^ for the rigged configurations in Example \3.2[ For 
convenience the vacancy p^'^^ is shown to the left of each block {a,j). The rightmost cell 
in the singular strings to be shorten are indicated by x . 

x]l oDo 











1 10 1 


4 




X 


4 





4 






3 
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(1") 
4 

(lio) 1 

(19) 2 
6 

(18) 2 



1 








3 




X 


3 



(1') l |x|x|x 
(1^) 



|0 




x]o oDo 



1 oDo 

1 O0O 

1 

1 



Thus the image is the highest path 11112221322433 G Bf^^. The algorithm of (p 
proceeds backward. 

Remark 3.6 Let V{L,\) = {p E Bf^\wtp = A} be the set of all the weight A paths 
including non highest paths. It is known that the algorithms for (p and (f)~^ actually 

work in a wider setting so that l{3.16\) is generalized to (f)(V{L, \)) ^ 'PiL, A). The set 

(j)(V{L, A)) of extended rigged configurations is characterized by l{3.1^) with a non-trivial 
lower bound ITW. 



3.3. Inverse scattering method 

Let p = hi ® ■ ■ ■ ® E Bf^ be a state of the BBS satisfying the boundary condition 
f l2.4ip . Suppose that p is highest and of weight A, i.e. p E V+{L,\). Then the state 
Ti{p) after the time evolution also belongs to V^{L,X). Thus, via the KKR bijection 
( ]3.16p . Ti on BBS states induces a time evolution of rigged configurations. The following 
theorem presents its explicit form. 

Theorem 3.7 [50i Prop. 2.6] For the subset of paths that undergo time evolutions 
without boundary effects, the commutative diagram 



RC(L,A) 



T, 



T, 

holds with the following time evolution Ti on rigged configurations: 



(3.17) 



RC(L,A) 



Ti : (/i, J)l h-> (/i, J')l, 

J'=(J'(1),J(2),...,JW] 



^!i^ = ^:^ + niin(/,j). 



(3.18) 
(3.19) 
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Namely, the KKR bijection linearizes the dynamics. More concretely, we find 

(/i^^), . . . , /i^")), ( J(2), . . . , are conserved (action variable), (3.20) 

Jj^^ changes linearly (angle variables). (3.21) 

Let us write fl3.19p as J' = J + hi, where hi = (5^1 min(/, j))^ plays the role of the 
(Zth) velocity vector. The commutative diagram fl3.17p provides a solution of the initial 
value problem in BBS. For a state p, it is given as T['{p) = (f)~^ o Tf^ o (j){p), where 
the T/^ in the RHS is just to change the rigging as J J + Nhi. The variety of time 
evolutions Ti, T2, . . . is refiected in the velocity vectors hi,h2, . . .. 

Example 3.8 The time evolution of the rigged configurations under Too corresponding 
to Example \2.11\ 

/iW ^(1) ^(3) 



1 



58^1 









□ 4 + 4t 




1 








10 + 3t 











15 + 2t 



1 no 



In Example 13. 8^ one notices that /i^^^ = (4, 3, 2) gives the list of amplitudes of 
solitons. This fact holds in general, which is a manifestation of the Bethe ansatz 
structure in BBS: 

/i*^^-* = list of amplitudes of solitons. (3.22) 

We call it soliton/ string correspondence. In fact, yU.^^^ is related to the earlier introduced 
conserved quantity Ei f l2.44p as 

El = number of cells in the left / columns of (3.23) 

There are still more conserved quantities in f l3.20p than They are responsible for 
the internal labels of colliding solitonfl. 

The inverse scattering scheme explained so far is naturally extended to not 
necessarily highest states by Remark 13.61 as long as the boundary effect is absent. For 
n = 1, the solution of the initial value problem in the same spirit has also been obtained 
in [68]. 

It was an essential insight of the quantum inverse scattering method [66] that 
Bethe ansatz can be viewed as a quantization of the classical inverse scattering method 
[H [19]. It is gratifying to realize that the combinatorial version of the Bethe ansatz 
here provides the inverse scattering scheme of the BBS which is a crystalline quantum 
integrable system. In this respect, the KKR maps (j) and are the direct and inverse 
scattering transforms and the rigged configurations play the role of scattering data [50] . 

+ The data (|3.20p is regarded as a rigged configuration for s[„ (instead of s(n+i) and tlie solitons are 
determined as tlie image of it under the KKR map (p~^. See [50l [54] for detail 
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4. Ultradiscretization — min-plus algebra 
4.I. Tropicalization and ultradiscretization 

Define T = M U {00} where 00 is the infinity which satisfies a < 00 and 00 + a = cxd for 
any a G M. The algebra (T, ©, ©) is called the min-plus algebra (or the tropical semifield) 
[63] , where the two operations "©" and "0" in T are respectively called tropical addition 
and tropical multiplication defined by 

a © 6 := min(a, b), aQb := a + b. 

The additive identity is 00, and the multiplicative identity is 0, i.e. 

a©cx) = a, a0O = a 

hold for any a G T. We have the inverse of © as a © (—a) = 0, but not the inverse of 
©. In the following we also write (T, min, +) for (T, ©, ©). 

We are to introduce a limiting procedure called the tropicalization, which links 
the subtraction-free algebra (M>o,+, x) to the min-plus algebra. We define a map 
Log. : M>o — M with an infinitesimal parameter e > by 

Log^ : a t-> —eloga. (4.1) 
For a > 0, define A G M by a = e~~. Then we have ljOg^{a) = A. Moreover, for a,b > 

A B 

define A, B E hy a = e~~ and b = e~~ . Then we have 

A B 

Log^(a + 6) = — e log(e"~ + e~~), Logg(a x 6) = A + 5. 

In the limit e — > 0, Logg(a + b) becomes min(A, 5). In this manner, the algebra 
(M>o, +, x) reduces to the min-plus algebra, and the procedure lim£_5.oLog£ with the 
transformation as a = e~~ is called the tropicalization. 

Through the tropicalization, subtraction-free rational equations on M>o reduce to 
piecewise-linear equations on M described by min-plus algebra, which is summarized as 
follows: for A, fi, C G M set 

_A _B _C 

a = e ^ , b = e <^ , c = e ^ 
and take the limit 5 — of the image Log^ of the equations 

(i) a + 6 = c, (ii) ab = c, (iii) = c. 

Then we obtain 

(i) mm{A, B) = C, (ii) A + B = C, (iii) A-B = C. 

We remark that the distributive law of the algebra (M>0) x)? o(^ + c) = a6 + ac, 
reduces to that of the min-plus algebra, A + mm{B, C) = m.m{A -\- B, A -\- C). 
Let us show an example. 
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Example 4.1 The discrete Lotka- Volterra equation for the variables {vj^ \ {j,m) G Z^} 
is given by 



V 



m+l 
j 



vf 1 + 5vf^^^ ' 



(4.2) 



where 6 is a positive parameter. We restrict vj^ G M>o, and take transformations 6 = e e 
and v"^ = e ~ . Then the tropicalization of (14. 2 p is calculated as 



J 

VJ^+^ - VJ" = -lim£(^log(l + e^^) - log(l + e 

= min(0, VJl^ + 1) - min(0, V;^' + 1). 



(4.3) 



By construction, the tropicalization of a discrete equation is defined on M, i.e. the 
dependent variables of the tropicalization are in M. When the tropicalization is defined 
on Z, we call it the ultradiscretization of the discrete equation. In the above example, 
(14. 3 p allows the ultradiscretization, since V™^^ is determined as an integer if VJ^, V"^i 



and V^^l"^ are integers. 



Remark 4.2 The original Lotka- Volterra equation vj = Vj{vj^i—Vj_i) is the continuous 
limit (5 — !■ o/ (14. 2 p with = Vj{—6m). Here 6 is a unit of the discrete time and vj is 
a derivation of Vj = Vj{t) by the time t. 



4-2. Evolution equations of BBS 

The original BBS in §1.11 corresponds to the time evolution Too in the formalism of §2.3[ 
which is the only case that admits the algorithms (albeit non-local) without carrier. 
One can set up two kinds of evolution equations for it: 

(i) the equation for the number of balls in the k-th box at time t [79] (the spatial 
description) , 

(ii) the equation for the number of balls in the j-th soliton (from the left) and the 
number W| of empty boxes between the j-th and the j + 1-th solitons at time t [77] 
(the soliton description). 

These descriptions are respectively related to the ultradiscretization of famous integrable 
difference equations, the discrete Lotka- Volterra equation ( §4.2.ip and the discrete Toda 
lattice equation ( §4.2.2p . 

4-2.1. Lotka- Volterra equation and infinite BBS. Let ul. G {0, 1} be the number of balls 
in the k-th box at time t. The evolution equation for u| is described by a piecewise-linear 
equation [79] : 

fc-i 

4+^=min(l-4, 5^(^-«n). (4.4) 



j=-oo 
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This equation has a piecewise-hnear version of the bihnear form in the following sense: 
assume that the variables {p^ | /c, t G Z} satisfy 

+ pi~' = ^^{pUi + pI pi+\ + pi''' - 1) • (4.5) 

Then the variables {m^ | /c, t G Z} defined by 

4 = pi + Pt\ - Pk' - Pi-i (4.6) 

satisfy (Oil . 

On the other hand, the discrete Lotka-Volterra equation (14. 2 p has a bilinear form: 

(1 + ^K+ir^^' = K+2rJU + ^T+i^T^ (4-7) 

i.e. if the variables {rj" | j, m G Z} satisfy the bilinear difference equation (14. 7p . then 
the variables {v^ | j, rra G Z} defined by 

m+l m 

satisfy fOj) . 

Proposition 4.3 [TS] i??. (14. 5 p zs i/ie ultradiscretization of the bilinear form (14.71) wi/i 
i/ie transformations 5 = e"^ anc? o"^ = under a coordinate transformation (j\ := T^_f 



Proof. It is obvious that (14. 5 p can be defined on Z. Via the coordinate transformation, 
(14. 7p becomes 

By applying the tropicalization with the transformation, we have 

- lime log(l + e"^) - p*" - p*+ ^ = - lime log(e ^ + e i ) , 
which yields (14. 5p . Here we use 

1 A. B 

limelog(l + e~") = 0, lim£log(e~ + e~) = max(yl, B). 

□ 

Remark 4.4 At (14. 4p . t(;e can regard v\_^ := X]j=-oo(^ "''^j''''^) '^■^ "^^^ number of balls 
in the carrier", which is identified with v^-i at (I2.39p . Then (14. 4p rewritten as 



u 



= min(l - 4, 4 = 4 + ^l-i " (4-10) 



One sees that these correspond to the description of BBS with the combinatorial R of 
crystal in ^2.2[ In fact, for n = 1, (I2.3ip simply reads 

Xi-Xi = -jji + yi = min(xi+i, yi) - min(xi, yi+i). (4.11) 

Thus by setting y = (1 — u\, u\), y = (1 — m^'^'^, m^^"^) G Bi and x = (9 — vl_i, vl_i), x = 
{9 - vl, vl) G Be, (I2.3ip reduces to (I4.10p in the limit 9 ^ oo. 
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The s[„+i BBS also has the bihnear form as fl4.5p . Given a state at time 
t = 0asp=---®x°(8) x1^^ ® ■ ■ ■ with x° G Bi, we consider its time evolution 
^(X)(p) = ■ ■ ■ ® ® xl^]^ ® ■ ■ ■ for t > 0. Here G -Bi specifies the local state of the 
k-th box at time t. According to f l2.17p . we express it as x^ = {x^ ^, . . . ,xl. n+i)- Define 
pliEZ{kEZ,i = 0,l,...,n+l,teZ>o)hj 



Pk,i — ^ (^j,2 + ^i,3 + ■ ■ ■ + X*j,i) + ^ ^ (2^j,2 + ^i,3 + ■ ■ ■ + ^j>+l) 
j=—oo t'>t+lj=—oo 



(4.12) 



Pfc,0 ~ Pk,n+1 

This counts the number of balls in the SW quadrant of the time evolution profile as in 
Example I2.11[ The variables {k, t) specify the position of the top right corner of the 
quadrant. The first term in (14.121) means that only those balls with color < i are counted 
on the top row of the quadrant. The quantities p^j are finite due to the boundary 
condition and the BBS time evolution rule. To see this concretely, note that each state 
at time t has a finite number of balls, therefore we have x* ^ = (j < /cq, i = 2, . . . , n + 1) 
for some feg ^ k. Then, the time evolution rule (Proposition l2.6p implies that the nonzero 
contribution to the double infinite sum in (14.121) actually comes from the finite region 
depicted as 



. . . t+2 

•^ko+2 -^k 



(4.13) 



^k 



t+k—ko 



By the definition, we have p^i^^+i — Pki 

4,i = pU - pLl,^ - PKi~i + 2 = 1, . . . , n + 1. (4.14) 
Proposition 4.5 [SH Prop. 4.2] The following relation holds: 

Pif-i + pLi,i = max(p*+/ + + pU "1) i = 2,...,n + l. (4.15) 

We note that the similar fact is studied in [21], IV]. When n = 1, we recover (14. 5p and 
(14. 6 p via Pk = Pk 2^ "^k — ^k 2- The variables p^. ^ will play an important role to solve the 
BBS in §4.4. 

4-2.2. Toda lattice and infinite BBS. Consider a state of the BBS with solitons, 
and let Q*- be the number of balls in the j-th soliton and Wj be the number of empty 
boxes between the j-th soliton and the j + 1-th soliton at time t as follows: 

Q\ Wl Q\ Wl VK^_i Q*v 

....1111111 52?25lL^£!^Tr^ IT^^J^iiiii 
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We have positive integers Q'j for j = 1, . . . , and Wj for j = 1, . . . , — 1, and have 
= Wj^ = oo. 

The evolution equations for and Wj are written as [77] 

j i-i 

Q';-'=mm{Y^Qi-Y,Qi+\ W^) j = l,...,N, (4.16) 

fc=i fc=i 

= g*+i + w^- g;+i j = i, . . . , iv - 1. (4.i7) 

Example 4.6 Lei «s consider the case of N = 3. In the following, the evolution of a 
3-soliton state at the left is written in terms of {Q\, W^, Q\, W2, Ql) at the right: 

t=0 . . .2222. . .222. . .2 (4,3,3,3,1) 

t=l 222... 222. 22 (3,3,3,1,2) 

t=2 222. . .2. .2222 (3,3,1,2,4) 

t=3 222.2 2222 (3,1,1,5,4) 

t=4 2.222 2222 (1,1,3,6,4) 

t=5 2... 222 2222... (1,3,3,7,4) 

One sees that the variables {Q{,W^,Q2,W2,Ql) satisfy i\A.16\i and (14.171) . 
On the other hand, the discrete Toda lattice equation is given by 

= «5 + u>5 - w'*\ (4.18) 

= (4-19) 

for j, t G Z. Now we only consider (I4.18P for j = 1, . . . , N and (I4.19P for j = 1, . . . , — 1 
with the boundary condition Wq = w]^ = 0. This is what is called the discrete Toda 
molecule equation. 

Proposition 4.7 [TTJ §3] Eqs. (14.161) and (14.171) are the ultradiscretization of the 
discrete Toda molecule equation with g* = e ~ and w^, = e ~ . 



Proof. By using (I4.19P iteratively, (I4.18P becomes subtraction-free: 



We apply the tropicalization and obtain the claim. Note that the boundary condition 
for Wq and is consistent with that for Wq and W^. It is clear that ( ]4.16p and ( 14.17^ 
are defined on Z. □ 

Remark 4.8 The original Toda lattice equation xj = e^-'+^~^-'' — e^-'"^^"^ is the 
continuous limit S of (I4.18P and (14.190 with Wj = S'^e^^+^~^^ and g* = 1 + Sxi, 
in the same manner as the Lotka-Volterra equation at Remark \4.2\ 
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Remark 4.9 The description fl4.16p . fl4.17p can be generalized to the infinite BBS of 



type s[.„+i J 77, §5y. Again we consider a state with N solitons, where we regard a 
non-increasing sequence of 2,3, ... ,n -\- 1 as a soliton ( ^2.3.^) . Let Q^j^ be the number 
of i-balls in the j-th soliton and Wj be the number of empty boxes between the j-th 
soliton and the j + 1-th soliton at time t. Then we have non-negative integers Qjj 



(j = 1,...,N, i = 2,...,n + l) andWj (j =0,...,N) satisfying ^"J2 Qli > 
Wq = = oo. The other Wj can be zero only if the color of the rightmost ball in the 
j-th soliton is strictly smaller than that of the leftmost ball in the j + 1-th soliton. 

We define Wj for t ^'Ljn, and regard Wj with t ^ Z as the intermediate states. 
The evolution equations are written as 

j i-i , 

^ V j.\ n-\-l — i 

Q5? = min(E^M-EC' K^) J = ^,---,N, (4.21) 

k=l k=l 

I I n-\-2 — i I I n + 1 — z 

^/ " =QUi, + W^^ " -Q5+^ i = l,...,iV-l, (4.22) 

where we run these equations from i = n-\-ltoi = 2. These piecewise-linear 
equations correspond to the ultradiscretization of the generalized (or hungry) Toda 
molecule equation. 

Remark 4.10 In this description only the information of relative coordinates of solitons 
survive, and the information of the absolute coordinates are lost. However, it is sufficient 
to study the basic features of BBS such as the soliton scattering and the conserved 
quantities. See fTT^ for the detail. 

4.3. Birational R and Geometric crystal 

The purpose of this subsection is to introduce birational R and geometric crystal for 
sin+i [IHl [52]. Besides their conceptual importance, they are useful to describe local 
evolution rules of discrete integrable systems related to BBS. They consist of birational 
maps and many other relations between certain sets of variables. The combinatorial R 
and the crystal for sl^+i in §2.21 are obtained from them by ultradiscretization. 

4.3.1. Birational R. Let B = {x = {xi, . . . ,Xn+i)} C (C^)"'^"^ be a set of variables. 
The birational R (introduced under the name of tropical R in [IHl IS2]) for is the 
birational map R: BxB^BxB specified by R{x, y) = {y, x) in which 

Pi-i{x,y) _ Pi{x,y) 



Pi{x,y) ' ^''Pi-i{x,yy 

n+l /n.+l k \ (4.23) 

k=l \j=k j=l 

where all the indices are considered to be in l^n+i- It satisfies the inversion relation 
R? = id on B X B and the Yang-Baxter equation 

R1R2R1 = R2R1R2, (4.24) 
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on B X B X B, where Ri{x,y,z) = {R{x,y),z) and R2{x,y,z) = {x, R{y, z)). A proof 
will be given later (Proposition I4.18|) . 

The birational R is characterized as the unique solution to a version of discrete 
Toda lattice equation. 

Proposition 4.11 [HH Th.2.2] Given {x,y), the birational R is the unique solution to 
the equations on {y,x): 

1111 , , 

XiVi = yiXi, \ = — + — , (4.25) 

Xi Ui+l Hi XiJ^i 

with an extra constraint nr=/(^j/^i) — nr=/(l/j/^i) — ^■ 

Proof. We prove that the ijj^x) given by (I4.23P satisfies (I4.25p . The uniqueness of 
the solution will be discussed later (Proposition I4.16|) . The former equation is clearly 
satisfied by fl4.23p . Let us check the latter one. It is equivalent to Pj+i (x, y) /xj+i + 
Pi_i(x, y)/yi = Pi{x, y)/xi + Pi{x, y)/yi+i which is verified as 

n / n k n k 

P,+i(x,?/)/xi+i - Pi{x,y)/xi = iYlxi+j+iYlyi+j+i - JJ x^+j JJ y^+j 

fc=l \j=k j=l j=k j=l 

n+1 /n+1 k n+1 k 

= ( n n y^+^ ~ n ^^+^-^ n y^+^-^ 

k=2 \j=k j=2 j=k j=2 

= Pi{x,y)/yi+i - Pi^i{x,y)/yi. 

□ 

In order to relate the birational R to the combinatorial R, we introduce the max- 
plus version of the tropicalization. It is a slight modification of the tropicalization in 
§4.11 For a > 0, define A EMhj a = e^. Then we have — Logg(a) = A. Moreover, for 

A B 

a,b > define A, 5 G M by a = eT and b = e~ . Then we have 

A B 

— Log^(a + b) = e log(e~ + e~), — Log£(a xb) = A + B. 

In the limit e — 0, — Log£(a + b) becomes max(y4, P). In this manner, the algebra 
(M>05 +5 x) reduces to the "max-plus" algebra, and the procedure — lim£_5.o Log^ with 

A 

the transformation as a = e~ is also called the tropicahzation. As in §4.11 it 
is called ultradiscretization when defined on Z. We note that this version of the 
ultradiscretization of (1423]) is ( 123T]) . and that of (14:25|) is (12:32|) . when we take 
B = {x = (xi, . . . G M>ofor alH}. 

4-3.2. Geometric crystal. A representation theoretical background for the birational R 
is provided by the geometric crystals [7] and their natural extrapolation into the aflfine 
setting [ini §1]. We explain this notion for sln+i- 

To give an overview of the basic idea, first we show a few relations in the case of 
up to 2- fold tensor products. Let us begin with the crystal. As a result of fl2.20p - (l2.25p 
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and by interpreting /j = e^^, the action of the Kashiwara operator with a parameter 
c G Zi is given, unless they vanish, by 

e- (x) = (..., Xi-i,Xi + c, Xi+i - c, Xi+2, ■■■), (4.26) 

e^(x®y) = er(x)®gf(y), 

ci = max(xi + c, yi+i) - max(xi, yi+i), (4.27) 
C2 = max(xi, ?/i+i) - max(xi, yi+i - c). 

In the geometric crystal, one still has the coordinates x = (xi, . . . ,Xn+i) G B and the 
corresponding structure looks as (c G C^) 

e- (x) = (..., Xi_i, cXi, c~^Xi+i,Xi+2, ■■■), (4.28) 

et{x,y) = ieTix),er{y)), 

cXi + yi+i Xi + yi+i (4.29) 

Xi + yi+i Xi + c ^yi+i 

We call el the geometric Kashiwara operator. Note that the ci, C2 in fl4.27p are piecewise- 
linear and obtained from f l4.29p by the ultradiscretization, i.e. replacing +, x,/ with 
max, +, — , respectively. 

Now we define the geometric crystal for s[.„+i in more general setting. In what 
follows, let c G be a parameter which takes generic values, el be a rational 
transformation on a variable set V C (C^)^ where G Z>o, and ei,7i are rational 
functions on V. 

Definition 4.12 A geometric crystal for sln+i is a family {V, £i,7i,e^} which satisfies 
the following relations. For any x E V, c,c' E , and z, j G / = {0, 1, . . . , n}, 

(i) e'ref{x) = ef{x), e}{x) = x, 

(ii) ei{e^i{x)) = c-^ei{x), 

(Hi) 7i(e^(x)) = c^7i(x) {i = j), = c~^7i(x) {i - j = ±1), = 7i(x) {otherwise), 

(iv) e^e^'(x) = e/e^(x) if i - j ^ ±1, 

(v) e^Cj'^'ef (x) = e'^j ef' e'^j{x) if i — j = ±1. 

Here i = j means i — j E {n + 1)Z. 

In what follows, we introduce the function ipi by 7j = ^i/Si. 

Example 4.13 For x E V = B, define el by ^4-^^ ^'^'^ ei{x) = Xj+i, (pi{x) = x,. 

Example 4.14 For {x,y) E V = B x B, define ef by ^4-^^ '^f^d let ei{x,y) = 
Xi+i(l + ?/i+i/xi), ^i{x,y) = yi{l + Xi/yi+i). 
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Example 4.15 For x = (x^^\ . . . , x'^^'^) G V = , define et, ^pi, el by 

eAx) = ^ r ^ 4.30 

n-=-iV.(x(^)) 

<P^i^) = ^ ^> (4-31) 

e^(a.) = (er(a;W),...,e^^(a;W)), 



with ci 



iJere 6'(s) = 1 if s is true and = otherwise. The Ei, ipi in the right hand sides are those 
defined in Example \4-13 . 

When c, c' G M>o, all the above relations in the geometric crystal for s[„+i reduces to 
the corresponding relations in the crystal for s[„+i in §2.2.11 via the max-plus version 
of the ultradiscretization. 



4-3.3. Matrix realization. There is a matrix realization of the geometric crystal for 
sln+i, where each element x G i3 is associated with the matrix 



M{x,C) 



■1 X'r, 



-1 



(4.33) 



1 ^n+lj 



involving the spectral parameter (. The structure of s[„+i geometric crystal is realized 
as simple matrix operations. The action of the geometric Kashiwara operator is induced 
by a multiplication of (product of) M with certain unipotent matrices. For simplicity 
we assume i 7^ in what follows. Let Gi{a) = E + aEi^i^i where E is the identity 
matrix. Then we have 



Gi 



c- 1 



EiiX 



M{x,C)Gi 



Lpi[X) 



M{e^{x),C), 



(4.34) 



for e^(x) G i3 in (14. 28 p . In the same way the action of on {x,y) G in fl4.29p is 
represented by 



G,, 



c- 1 

£i{x,y) 



M{xX)M{yX)G^ 



- 1 



^i{x,y) 



M{ef{x)X)M{eT{y)X). (4.35) 



By using the formulas in Example I4.15[ one can also define the action of for 
general multiple product case. Through these examples, we observe that the product 
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of matrices M {x^^\ ■ ■ ■ M {x^^\ corresponds to the product of geometric crystals 

It is easy to see that equation f l4.25p is equivalent to the matrix equation 



M{x, C)M{y, = Miy, C)M(5, C). (4.36) 

Due to the presence of the spectral parameter (, its non-trivial solution is unique. 
This characterizes the birational R as the unique intertwiner (i.e. the operator that 
interchanges the order of product) of the geometric crystals. For x = (xi, . . . , Xn+i) € B 
we set £{x) = xi ■ ■ -Xn+i- The uniqueness of the solution to ( 14.36^ can be verified as a 
consequence of the following: 

Proposition 4.16 Suppose M{x,C)M{y,() = M{x'X)M{y\Q for i{x) = i{x') ^ 
£{y) = £{y') . Then x = x',y = y' . 

Proof. We define M{xX) = (1 ~ C^{.x))^{.^X) to avoid a singularity of M(x, C) at 
( = £{x)-^. Now the relation M(x, C)M(y, C) = M(x', C)M(?/', C) is supposed to be 
satisfied under the condition £{x) = £{x') ^ £{y) = £{y'). It is easy to see that the 
matrix elements of M(x, Q are given by 



C \\[k=i^kj \ \[k=j Xk) for z < J. 



Note that the rank of the matrix M{x, () reduces to one when ( = i{x) ^ or more 
precisely M{x, i{x)~^)ij = i{x)~^ (^111=1 ^fcj (^Ilfeij ^fc j for any Thus the relation 



M{x,C)M{yX) = M{x',C)M{y',C) at C = ^(x)"^ yields the condition (Jll^^Xk 

a ^ni=i ^'k^ for any i and with some constant a, forcing Xi = x[ for any i. In the same 
way we obtain y = y' hj taking ( = i{y)~^ in the relation. □ 
Now we show that the birational R is the intertwiner of the geometric crystals. 

Proposition 4.17 Re^ = e^R. 

Proof Let R{x,y) = (y,x), e'rR{x,y) = ((y)', (x)'), e^(x,?/) = {x',y') and Re1{x,y) = 
{y',x'). Then we have 

Af(y', C)M(?, C) = M(x', OM(y\ Q 

\'^i\-^^y)/ \H^i\-^jy) 



^ei{y,x)J ' ' \ipi{y,x) 

= M{{yyX)M{{x)'X)- 

Here we used ei{y,x) = ei{x,y) and (pi{y,x) = ipi{x,y) which are verified by fl4.25p . By 
Proposition 14.161 we have y' = {yy,x' = (x)'. □ 
We show that the birational R satisfies the Yang-Baxter equation. 
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Proposition 4.18 P", Th.2.2] The birational R satisfies 

Proof. Let RiR2Ri{x,y, z) = {z',y',x') and R2RiR2{x,y, z) = {z",y",x"). By 
Proposition 14.111 and since f l4.25p is equivalent to f l4.36p we have 

M{z',C)M{y',C)M{x',0 = M{x,OM{y,C)M{z,C) = M{z" ,C)M{y" ,C)M{x" ,(). 

(4.38) 

By an obvious extension of Proposition 14.161 this leads to x' = x", y' = y", z' = z". □ 

4.3.4- Bilinearization. The birational R is equivalent to a system of bilinear difference 
equations of Hirota type [52]. To see this, introduce the functions r/ (1 < J < 4,i G 
Zi„+i) and the parameters Aj, Kj, and make the change of variables 

xj^ = Xi^rf I drf , yj'^ = Kidrf / dr} , 
\ V 4 -1 4 1 (4-39) 

with 5t^ = T'l ItI_y. In order to memorize the relations fl4.39p it is useful to draw the 
following vertex diagram and regard the tau functions as residing in the quadrants. 



X- 



; 








; 





(4.40) 



Then the former relation in f l4.25p is automatically satisfied and the latter is translated 
into 

XiTtirt - ^^ir■Tt_^ = arlrl, (4.41) 

for any nonzero parameter a independent of i. The birational map R : {x,y) ^ iy,x) 
is induced by an automorphism ^ t^, \i ^ Hi, a ^ —a of (14.410 . Eq. (I4.4ip is a 
version of so-called Hirota-Miwa (non-autonomous discrete KP) equation. 

4.4- General solution 

Recall that the KKR map ( §3.2p transforms a rigged configuration into a highest 
path. It turns out that its image allows an explicit formula in terms of ultradiscrete tau 
functions. In view of the remarks after (13.230 . this yields the general solution of the 
BBS time evolution equation (I4.15P corresponding to an arbitrary initial condition. 

To formulate the ultradiscrete tau function, it is convenient to regard a rigged 
configuration ( §3.2p as a multiset i.e. a set with multiplicity of each element taken into 
account 

S = {(fli, Ji) e {1, 2, . . . , n} X Z>i X Z>o M = 1, 2, . . . , A^}, (4.42) 
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where > is arbitrary and each triplet s = {a,l,r) signifies a string having color a, 
length I and rigging r. This fact will be denoted by cl(s) = a, lg(s) = / and rg(s) = r0- 
S" is a rigged configuration if rg(s) < p{g(^')^'' is satisfied for all s G S", where p^"^ is defined 
in ([31D. 

For a rigged configuration S (14.421) . let T C 5 be a (possibly empty) subset of S. 
We allow the fact that T is no longer a rigged configuration in general. Introduce the 
piecewise-linear functions Ck^a{T) {0 < k < L and l<a<n + l)by 

<T) = IY1 Cci(.),ciwmin(lg(.),lg(t)) + J]rg(.), (4.43) 

s,teT seT 

Ck,a{T) = c{T) + Yl Us)-k Y 1' (4.44) 

s£T,c\{s)=a sGT,c1(s)=1 

where Ca,b = '^Sa,b — <^|a-6|,i is an element of Cartan matrix {Ca,b)i<a,b<n of s[„+i. (L will 
be the length of the corresponding path.) By the definition, the second term in Ck^a{T) 
is when a = n+1. Obviously we have c(0) = Cfc^a(0) = 0. The quantity c{S) is known 
as the cocharge of the rigged configuration S 03]. The ultradiscrete tau function is a 
Z>o-valued piecewise-linear function r^^a = 'Tk,a{S) on S defined by (0 < /c < L) 

rfc,a = -min(cfc,a(T)) (l<a<n + l), 

(4.45) 

Tfc,0 = T-fc,n+l - k. 



Example 4.19 To,„+i = — minTcs(c(T)) for S i4-42^ with = 1,2,3 is given by 

7-o,n+i = -min(0,^i), 

i"o,n+i = - min(0, ^i, ^2, 6 + 6 + ^1,2), 

ro,n+i = - min(0, 6, 6, 6, + 6 + ^1,2, ^1 + ^3 + v4i,3,6 + 6 + ^2,3, 
6+6 + ^3 + ^1,2 + ^1,3 + ^3), 

where we have used the shorthand 6 = li + o't-'^ ^i,*' = Ca^^a-, niin(/j, /j/). 

In general, the minimum (I4.45P for S (14.421) extends over 2^ candidates and reminds 
us of the structure of tau functions in the theory of solitons pD]. In fact, (14.451) can be 
deduced from the tau functions in the discrete KP hierarchy by ultradiscretization with 
an elaborate turning of parameters between KP solitons and rigged configurations [SU 
sec. 5]. 

In §3.31 we have seen that rigged configurations undergo linear time evolution (I3.19p . 
In the present notation, it is rephrased as 

S = {(a„ k, J,)} A Ti{S) := {(a,, J, + min(/, /,))}• (4.46) 

* Colors 1,2, ... ,n of strings in rigged configurations should not be confused with the colors of balls 
in BBS. 
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Theorem 4.20 Let 61 (>>>■■■ ® 6l = 4'~^{S) be the image (highest path) of a rigged 
configuration S under the KKR map (f>^^ . 

(i) [SU Th. 2.1] hk = . . .,Xk,n+i) e Bi (12.171) is expressed as 

Xk,a = Tk,a - Tk~l,a " Tfc,a-1 + T-fc_l,a-l- (4.47) 

(ii) [Ml Prop. 5.1] Denote by Tk,a the ultradiscrete tau function associated with T^oiS) 
defined by (Thus r^^i = T^^n+i-) Then the following ultradiscrete Hirota-Miwa 
equation is satisfied. 

Tk,a-1 + Tk^l.a = m.ax(Tk,d + rk^l,a~l, rfc_l,a-l + Tk^a "1) (2 < fl < n + 1). (4.48) 

(Hi) [SI Th. 4.9] Define p = • ■ ■ (g) (g) xl^^ (g) ■ ■ ■ by = bk if 1 < k < L and x° = |T| 
otherwise. Let pi ^ be the number of balls specified from p as in ( [^.i^ ^. Then Tk,a = pt.a 
holds for l<k<L,l<a<n-\-l. 

Theorem 14. 20l is known to hold also for extended rigged configurations (Remark 13. 61) 
and non highest paths [Ml sec. 7]. In view of the inverse scattering method ( §3.3p . it 
provides the explicit piecewise-linear formula describing BBS under any time evolution. 

Remark 4.21 Let Na be the number of strings in a rigged configuration S having color 
a. The soliton/ string correspondence liS. 22) tells that S describes the Ni-soliton states 
of BBS. On the other hand, Tk,a{S) is an ultradiscretization of an N = Ni + ■ ■ ■ + Nn- 
soliton solution of the discrete KP equation \5^. For n > 1, the "extra" N2 + ■ ■ ■ + 
solitons in KP specify the internal labels of the BBS solitons. 

The cocharge mentioned under (I4.44p is related to the energy of a path, which 
involves the energy function H (I2.33P as a building block. See for example [EH [621 165] . 
In this context, ultradiscrete tau functions are combinatorial analogues of corner transfer 
matrices in solvable lattice models [6], and ( 14.47P is regarded as the formula for a "one- 
point function". 

These features and the insights gained in §3. 31 are summarized in the following table. 
One can compare the format of the solutions of BBS coming from the two basic tools 
in quantum integrable systems, Bethe ansatz and corner transfer matrices. 





Bethe roots 


Corner transfer matrix 


Combinatorial analogue 


rigged configuration 


energy in crystal 


Role in BBS 


action-angle variable 


tau function 


Dynamics 


linear 


bilinear 
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5. Periodic BBS 

5. 1 . Basic features 

In this section we restrict ourselves to s\2 case and consider the box-ball system with 
the periodic boundary condition, which we call periodic BBS for short [571 EHl ESI IM] . 
An example of the time evolutions of this system appeared in §1. Compared with the 
infinite system, there are many interesting features in periodic BBS which come from 
the finiteness of its phase space. For an attempt to generalize our formalism to the case 
of periodic sl„_|_i BBS, see [56] . 

Let us recall the formalism in §2.3l which is based on the crystal base theory. In the 
case of SI2, the vertex diagrams for combinatorial R fl2.37p on Bi ® Bi look like those in 
fl2.14p . Let L be the system size and M (<L/2) be the number of balls. Consider the 
diagram fl2.39p with not necessarily large L. Let 



VL = {h®---®hLe =\2\}< L/2}. (5.1 



To attain the periodic boundary condition we want to find v' G Bi such that v' = v. See 
fl2.15p for an example. In general v' is a function of v G Bi andp := 61®- ■ -^bi G 
Hence we can denote it by v' = v'{v,p). 

Proposition 5.1 [53, Proposition 2.1] For any p G Vl under the condition = 
2} < L/2, the solution v G Bi to the equation v'{v,p) = v is unique and is given by 



V = v'{ui,p), where = 1 ■■■ 1 G Bi. 

Let vi{p) = v'{ui,p). By setting v = v' = vi{p) in fl2.39p we define the time evolution 
operator Ti by the relation f l2.43p . That is, we have 

vi{p)®pc^Ti{p)®vi{p), (5.2) 

as elements of Bi ® (Bi)^^ ~ (i?i)®-^ ® Bi. We note that (15.20 is a periodic version of 
the Lax equation (12.451) . In particular Ti yields a cyclic shift by one unit cell to the 
right. The evolution by Too admits the description without carrier given in §1, which is 
also equivalent to the "arc rule" in |83j . 

Remark 5.2 In §5. il we restrict ourselves to the case M < L/2 for simplicity. However, 
our formalism of the periodic BBS based on the crystal base theory also enables one to 
treat the case M > L/2 fS^. 

The energy associated with T/ is defined by (I2.44p with vq = vi{p), where the values 
of the energy function are given hj H = for the bottom right diagram in (I2.14p and 
H = 1 otherwise. 

In what follows, we write for example [T] (S) [T] (S) simply as 1212. 
Example 5.3 The time evolutions of p = 222111211111 by Ti with / > 3,T2 and Ti: 
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t=0 222... 2 

t=l ...222.2 

t=2 2.222. 

t=3 22 2... 2 

t=4 ..222... 2... 

t=5 222.2.. 

t=6 2 2.22 



222. . .2 

. .222. .2 

222.2. . . 

22.22. 

2 2.22 

222 2. . 

. .222 2. 



222. . .2 

222. . .2. . . 
222. . .2. . 
.222. . .2. 
. .222. . .2 
. . .222. . .2 
222. . . 



Theorem 5.4 [571 Th.2.2] The commutativity TiTk{p) = TkTi{p), and the conservation 
of the energy Ei(Tk{p)) = Ei{p) hold. 



Proof. Let R{vj^{Ti{p)) ® vi{p)) = vi{p) ® Vk{Ti{p)). See the following diagram 



Vk{Ti{p)) 




VkiUp)). 



p 














Tlip) 














TkTiip) 



Mp) 

Vk{Ti{p)) 



Vk{Ti{p)) 



vi{p) 



p 














z 














w 



X 




vi{p) 
Vk{Ti{p)) 



where the Yang-Baxter relation is used to move the symbol "x" for the R from the left 
to the right. Consider what x, y, z and w should be. Since R{yk{Ti{p))®vi{p)) = vi{p) ® 
VkiTiij))) we have x = Vk(Ti{p)) and y = vi{p). Hence by Proposition 15.11 the equality 
VkiTiij))) = Vkip) holds. Thus z = T/^^p). Then again by Proposition 15.11 the equality 
vi{p) = viiTkip)) holds, and we have w = Ti{Tk{p)). Hence Tk{Ti{p)) = Ti{Tk{p)). For a 
proof of the conservation of the energy, see (STJ Th.2.2]. □ 



5.2. Linearization and general solution 

Here we construct action-angle variables of the periodic BBS, solve the initial value 
problem and present an explicit formula for 7V-soliton solutions in terms of tropical 
Riemann theta function. These results were firstly obtained in [511 [57] . 



5.2.1. Action variable. We are going to introduce the action variable of a state. It is 
equivalent to the list of amplitudes of solitons contained in a state, which is the conserved 
quantity. Recall that a state p = 6162 • • • &l [bj = 1, 2) is highest if the condition (13.41) is 



Integrable structure of box-ball systems 



49 



satisfied. It is elementary to show that any state in Vl can be made highest by a cychc 
shift. Namely, there exist a highest state G Vl and G Z such that p = Tf{p^). 
Given p, such a pair {d,p+) is not necessarily unique. Nevertheless one can show 

Proposition 5.5 [571 Proposition 3.3] (i) Let fi be the configuration of the rigged 
configuration (l){p+). (Namely, /i is the Young diagram denoted by /i^"^^ in ^3.2[ ) Then 
/i is independent of the not necessarily unique choice of {d,p+). 
(a) The energy Ei of p is related to /i via 



Ei{p) = number of cells in the left I columns of fi. 



(5.3) 



Due to (i), the Young diagram /i is uniquely determined from a state p. We denote 
it by /i*(p) and call it action variable of p. Due to (ii), it is a conserved quantity, namely, 
H^^iTiij))) = n^^ip) holds for any /. Let ruk = mk{p) be the number of length k rows in 
the Young diagram fi*{p). Then (15.31) is rephrased by the same formula as (12.471) : 



Ei{p) = y^min(/, k)mk. 



(5.4) 



k>l 



In the context of the KKR bijection, is the number of length k strings. On the 
other hand, rrik is the number of amplitude k solitons contained in p. This is another 
manifestation of the soliton/ string correspondence, which was observed earlier also in 
BBS on infinite lattice in (I3.22p . We introduce the isolevel set of states characterized 
by the action variable 

VM = {peVL\ Mp) = /i} (5.5) 



for any Young diagram fi such that < L/2. 

Example 5.6 Take p = 2211221112122111221 G Pig. 

cyclic shifts of highest paths as p = T^{p+) = Tf{j)'_^) 
1122111212211122122, p'^ = 1112122111221221122 and pi 
Their image by the KKR map (f) is given by 



It can be expressed as 
- Ti^ip'l), where p+ = 
= 1112212211221112122. 



P+ ^ 



4 













3 






1 




6 






2 





p'l 



They all lead to fi*{p) = (3,2,2,1,1). 
Example 5.7 The isolevel sets Vq{^) with \^\ 



3 are given by 



P6((l, 1,1)) = {121212,212121}, 

pg((3)) = {111222,211122,221112,222111, 122211, 112221}, 
Vfi{{2, 1)) = {121122, 212112, 221211, 122121, 112212, 211221, 
112122,211212,221121, 122112,212211, 121221}. 



(5.6) 
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5.2.2. Angle variable. Let us observe Example I5.6[ Recall that in the infinite system, 
the time evolution Ti is the uniform shift of all the riggings by 1. See f l3.19p . If we 
adopt the same feature also in the periodic BBS, the following identification should be 
made for the "periodic version" of the rigged configuration: 













3 






2 




10 




6 





12 



1 

= 3 
9 







|13 






16 






15 




21 




16 



(5.7) 

Here we have attached the vacancy pi on the left of the block of length i strings. The 
riggings are no longer bounded by it. The basic idea in constructing angle variables is 
to introduce an appropriate equivalence relation among such extended riggings. 

We proceed to the precise definition. Consider the iso level set Vl{ij) with /i = 
(i^'" ...i^'^). Here ii < ■■■ < ig are the length of the rows in fi and rrii. is the 
multiplicity of ij^. For instance /i = (3^2^1^) in (15.71) . We set 



X = {ii < ■■■ < ig}, pj = L - min(j, i)mi {j G Z>o) 



(5. 



where the latter is the vacancy pj^'^ (13. 9p with n = 1. Recall that a rigged configuration 
(/i, J)l is the data in which the vicinity of the block of length i strings looks as 



Pi 



J = (A 



z G X, 1 < a < m,-, 



< < ■ ■ ■ < Ji,m, < Pi 



(5.9) 



We extend the integer sequence Ji^a from 1 < a < mj to a G Z by imposing the quasi- 
periodicity as 

Ji,a+mi = Ji,a+Pi (« ^ Z). (5.10) 

The resulting sequence will be denoted by J = {Ji,a){ia)£ixz and called a quasi-periodic 
extension of the rigging J. (Indices will be suppressed as J = (Ji,a).) By the definition, 
J ranges over the set 



= JJ A(m,,pi), 

A(m,p) = {{Xa)aeZ I G Z, Ac, < Xa+l, Ac 

where L-dependence enter (15. lip via pi 



Xa+P (Va)}, 



(5.11) 
(5.12) 



jj The ii, . . . ,ig here and in ^6.31 will denote the amplitude of solitons. They should not be confused 
with the ones in §3] like p = ii ■■■ il or in p.4p and (|3.7p . 
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Now we introduce the equivalence relation on ^/^(/i). For G T, define (Xfc by 

CTk-JM-^JM; (Aa) ^ {Ji,a+5,„,+'2ram{i,k)). (5.13) 
Let A be the abelian multiplicative group generated by (Xj^ , . . . , (Tj^ . Define 

JM = JM/A (5.14) 

which is the set of equivalence classes of Jl{ij) under A. The image [J] G Jl{i^) 
of J G Jhilj) will also be written as J for simplicity unless emphasis is preferable. 
Elements of J^iifJ') are called angle variables. 

Angle variables are also depicted as (15. 9p . Actually, infinitely many such diagrams 
that are transformable by A all correspond to a single angle variable. For instance in 
(15. 7p . if the leftmost one is [J], then the middle and the rightmost ones are [o"2(J)] and 
[(Ti(t|(J)], respectively. 

We introduce the time evolution T; (/ > 1) on j7L(/i) by 

Ti : JM ^ JM; (Aa) ^ iJi,a + min(^, /)), (5.15) 

and denote its induced action on J^iifJ') also by Ti. Obviously, Ti is linear and 
commutative. In particular we use the abbreviation Tf{3) = J + (i for the uniform 
shift. Readers are highly recommended to check that (riigj^I"') (J) = J + G Jl{.Ij)i 
which implies that any angle variable is invariant under as it should. 

5.2.3. Linearization of time evolution. Let us assign an angle variable to each state 
in the isolevel set ViifJ')- Namely, we are going to construct a direct scattering map 
$ : Vl{ij) — Jl{,Ij)- We do this by suitably adapting the KKR map (p to the periodic 
setting. Let V^{fi) = {p G ^^(/i) \ p : highest} be the subset of Vl{ij) consisting of 
highest paths. We consider the following scheme: 

ZxP+(/i) ^ ^ JM ^^^g^ 

p I — )■ {d,p^) I — > 3 + d I — )• [J + c/]. 

First arrow: Pick any {d,p^) such that p = Tf{p^). Second arrow: Apply the KKR map 
(p{p+) = (/i, J)l and quasi-periodically extend the so obtained rigging J to J followed 
by a uniform shift by d. Third arrow: Take the image in J^iifJ') (identify by A). In order 
to make sense of the scheme (I5.16P as a definition of the map $, the non- uniqueness in 
the first arrow must be canceled in the identification in the third arrow. It was indeed 
the case in the example (15. 7p . Here comes the main result of this section. 

Theorem 5.8 [571 Th.3.11] $ is well-defined, bijective and makes the following diagram 
commutative. 



VM JM 



Ti (5.17) 

VM JM. 

Here Ti on the left and right sides are defined by ( f5. S\} and ( f5. iJj) . respectively. 
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The commutative diagram fl5.17p is the periodic version of fl3.17p . According to Theorem 
15. 8[ the nonhnear time evolution on Vl{h) is transformed to a straight motion on 
J^l{iJ')- This is a characteristic feature in finite dimensional integrable systems, where 
the dynamics is linearized, via what is called the eigenvector map, on Jacobian variety 
(or more precisely, abelian variety). In our setting here, the modified KKR map $ and 
the set of angle variables ^/^(/i) play the analogous roles to them. In §6.3 we will study 
this feature with tropical geometry, and see that Jl{ij) corresponds to the lattice points 
of some tropical abelian variety. 

Theorem 15.81 led to the first complete solution of the initial value problem of 
the periodic BBS. It is obtained by going along the commutative diagram fl5.17p as 

= o Tf^ o $. The variety of time evolutions Ti,T2,... is reflected in the 
corresponding velocity vectors in ( ]5.15p . 

Example 5.9 Let us derive a time evolution of the length 19 path p in Example \5.(k 

T^{p) = 1221112211211221122. (5.18) 



The angle variable of p has been obtained in ([5. 7| ). Using the leftmost representation, 
we find 



p I — y o 
9 













3 






2 




10 




6 





I > 



13 
12 



15 
11 



I — > 



|0 



The vacancies are exhibited only in the leftmost diagram. The rightmost diagram is 
a rigged configuration and corresponds to the highest path p' = 1121122112212211122. 
Therefore the image of the rightmost angle variable by is Tf{p') giving the RHS of 

(E23j. 

As this example indicates, to compute the inverse image $~^(J), one first finds a 
representative of the angle variable J of the form (rigged configuration) + e with some 
e e Z. Namely, one transforms J into C3"( J) = J'+e with an appropriate element a G ^ so 
that J' becomes the quasi-periodic extention of some rigged configuration (/x, J')l- Then 
one applies the KKR map to get a highest path as p'j^ = 0~^((/i, J')l)- Finally, 
the inverse image is obtained as $~^(J) = Tf (p'_,_). The fact that these procedures are 
always possible and the result is unique is guaranteed by Theorem 15. 8[ We note that 
the solution of the initial value problem based on the procedure called 10-elimination 
[58] is equivalent [H] to the preceding solution [57] explained here. 



5.2.4- N-soliton solution. Let us present an explicit formula of the path p = $^^(J) G 
Bf^ that corresponds to the given angle variable J G Jl{iA- This is a combinatorial 
analogue of the Jacobi inversion problem, and the result is indeed expressed in terms 
of tropical Riemann theta function. For simplicity we restrict ourselves to the case 
rrii = 1 for all i G X. (See [55] for the general case, where the tropical Riemann theta 
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function with rational characteristics is involved.) We retain the notation (15. Sp . where 
the latter reduces to pj = L — 2 J2iei ^^^U^ 0- The angle variable J = {Ji,a){ia)€Xxz can 
be presented simply as a ^f- dimensional vector J = {Ji)i^x with Jj = Jj^i since the other 
components are specified as Ji^a = Ji + {a — l)pi by the quasi-periodicity. It is easy to 
translate the identification by f l5.13p into this representation of J, leading to the simple 
description 

JM = Z^/FZ\ F = {5,,p, + 2min(z, j)),_^.gj. (5.19) 

The g X g integer matrix F is positive definite and has the origin in the study 
of Bethe equation at g = |18]. We introduce the tropical Riemann theta function 
(Definition iQll 

e(Z;F) = mm {n- (iFn+ Z)} (Z G M^) (5.20) 

in which the F is built in as the period matrix. We further introduce the (yf-dimensional 
vectors 

P = {Pi)iex, h = (min(i,/))^g^, (5.21) 
where the latter is the velocity corresponding to Ti fl5.15p . 

Theorem 5.10 [511 Th.3.3] The state ^~^{J) = 6162 ■ ■ - bi {bj = 1,2) corresponding to 
the angle variable J G J^ifJ^) is expressed as (0(Z) = 0(Z; F)) 



fefc = 1 - e (j - I - A;hi) + e (j - I - (A; - l)h, 

+ Q (3 - ^ - khi + h^) - e (3 - ^ - {k - i)hi 



2 — ■ -^y "V 2 



(5.22) 



The time evolution by Ti is attained just by replacing J with J + h;. Thanks to the 
quasi-periodicity (16. 6p . bk is invariant under the change J t— )■ J + FZ^. The invariance 
bk = bk+L is due to Fhi = Lhi. 

Remark 5.11 The matrix F f l5.19p will be related to the period matrix VL f l6.17p of the 

tropical spectral curve of the tropical periodic Toda lattice at Proposition \6.21[ 



5.3. Decomposition into Torus 

5.3.1. Introduction. Let us discuss the structure of the isolevel set of the periodic 
BBS. Recall the definition of the isolevel set 'PL(/i) in (15.50 . We use the notations 
H = {i^'^ ...i™'^), X, Pj in §5.2[ In contrast to the tropicalization of the periodic 
discrete Toda lattice in §6.21 where the isolevel set Tc given by (I6.14p is isomorphic to a 
real torus M^/fiZ^, the set Pl(/u) is a finite set. As we have shown in Theorem 15.81 the 
set 'PL(/i) is identified with the set Jl{,Ij) which was defined as a quotient set (I5.14p . 

The set Vl{h) can be regarded as a graph in the sense of graph theory. Let T be 
the abelian multiplicative group generated by Ti,T2, . . .. Then T acts on the isolevel 
set Vl{ij)- If one represents the elements of Vl{ij) by nodes and the actions of the time 
evolutions by arrows, then one has a colored oriented graph. The graph for ViifJ') has 
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usually several connected components. From the viewpoint of the group action we say 
that the isolevel set Vl{p-) decomposes into T-orbits. 

To begin with let us illustrate a few simple examples. First let us assume mj = 1 
for all j e I. Then we have = Z^/FZ^ (15.191) which is the set of all integer points 

on the torus MP / FX^ . We shall call IP jFV itself a (gf-dimensional) torus for short. 

Example 5.12 The isolevel set V^{{2)) is depicted as follows. 



,22111 



21112 



12211 




(5.23) 

Here the actions ofTi andT2 are represented by thin and thick arrows respectively. Note 
that P5((2)) ^ J75((2)) = Z/5Z. The velocity vectors Ii5.21\) are given by hi = (1) and 
h2 = (2) which reflects the relation T2 = (Ti)^ on this isolevel set. 

Now we consider the case with mj > 1 where the graph for Viif^) has indeed several 
connected components. Let 



S(p) = {p' E Vl{ij)\p' = gp for some g E T} C Vl{i^) 
be the T-orbit of p G Vl{ij)- 

Example 5.13 The graph for the isolevel set Vq{{1, 1)) is depicted as follows. 

y21211K 



(5.24) 



121 



12 



121211 



211121 112121 
^111212^ 



211211 



112112^121121 



(5.25) 



Here the actions of Ti are represented by arrows. It is not necessary to consider Ti>2 
because they effectively coincide with Ti in the present case. The graph has two connected 
components, as V^{{1, 1)) = S(212111) UE(211211). Each connected component has the 
structure of a one- dimensional torus, S(212111) ^ (Z/6Z) and S(211211) ^ (Z/3Z). 



5.3.2. Internal symmetry. The isolevel set in Example 15.131 has two connected 
components with different sizes. The difference reflects the internal symmetry. If the 
state has a larger internal symmetry, then it belongs to a smaller connected component. 
Let us briefly discuss this notion here. 

The internal symmetry of a BBS state is represented by an integer vector 7 = 
(7jj, . . . , 7j^) E (Z>o)^- We demonstrate how one can read off 7 from p. Given 
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p G ViifJ') there exist a highest path and G Z such that Tf{p^) = p. By the KKR 
map in §3.21 a rigged configuration is constructed from p^ as (f){p+) = {ii,J)l where 
J = {Ji,a)i£X,i<a<mi- We adopt the quasi-periodic extension for the angle variables 
fl5.10p . Now the integer & X fl5.8p ) is defined as the largest common divisor of 
and Pi such that 



J, 



i,a+- 



Ji,a ~l~ 



Pi_ 



[a G Z) 



(5.26) 



Neither the action of in (15.131) nor that of T; in (I5.15P changes the relation (I5.26p . 
Therefore Theorem 15.81 ensures that the internal symmetry of p is uniquely determined 
from the above procedure even if there is non-unique choice of p+. 



Example 5.14 Take p = 1212111222, p = 1211121222 G Vio in Example [XTTj They 
are already highest. By the KKR map we obtain the following rigged configurations. 



P ^ 4 






















P ^ 4 









2 












(5.27) 

In this example one has g = 2, X = {1,3}, hence 7 = (71,73). Note that the 
configuration implies = 1 which imposes 73 = 1. While Jia+i = Ji,a + | for p 
, no such symmetry exists for p. Hence one has 7 = (1, 1) for p, and 7 = (2, 1) for p. 

Example 5.15 Take p = 121122111212211222121111 with L = 24 wMch is already 
highest. By the KKR map we obtain the following rigged configuration. 



p I— )■ 



4 
6 

12 







(5.28) 



Here we have g = 3 and X = {1,2,3}, hence 7 = (71,72,73)- Observe that 
■/i^Q+l = Ji,a + y which implies 71 = 3. No such symmetry exists for J2,q, (and Js^a)- 
Hence one has 7 = (3, 1, 1). 



Let 



'PL,'y{fJ') = {p G 'PL(At)|the internal symmetry of p is 7}. 
Then we have ^^(/i) = U7'^l,7(/^)- 



(5.29) 



5.3.3. Connected component as torus. So far the states of periodic BBS are classified 

as 

Vl D VM ^ VlM ^ S(p). (5.30) 

Now we study the structure of a single connected component S(p). Let he a g x g 
matrix defined as 



F^ = F ■ diag(7-\...7-'), 
F = {6ijp^ + 2mm{i,j)mj)..^^. 



(5.31) 
(5.32) 
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This matrix is a generalization of F in f l5.19p . Then the structure of a connected 
component is stated as follows. 

Proposition 5.16 [TTl Th.2] Every connected component of the isolevel setViifJ') 
the structure of a g-dimensional torus "ZP jF^'ZP . The time evolution Ti is realized as the 
straight motion with a constant velocity vector hi = (min(z, on the torus. 

This result may be viewed as an ultradiscrete analogue of the classical Arnold- 
Liouville theorem [3]. 

Example 5.17 Take p = 1212111222, p = 1211121222 e Vio- They belong to two 
distinct connected components of the level set Pio((3, 1, 1)). The T-orbits and 
have the structure of two-dimensional tonZ^/F(i i)Z^ andl? / F(2,i)'^'^ respectively, 

where = \. ^ ) and -F(2,i) = ( o e ) • They are depicted as follows. 




(5.33) 



The nodes within and on the edges of the parallelograms represent the states of periodic 
BBS which belong to each of the connected components. Every pair of the parallel edges 
of each parallelogram should be identified. Thin and thick arrows represent the velocity 
vectors hi and \i2 respectively. 



5.3.4- Fundamental period. The time evolution T; is invertible and the isolevel set 
Vl{ij) is a finite set. Therefore every path p G Vl{ij) possesses the property T/^{p) = p 
for some integer > 1. We say any such integer a period of p. The minimum period 
is called the fundamental period of p and denoted by Mi{p). Every period is a multiple 
of the fundamental period. A formula for the fundamental period under any Ti was 
established in [57] . Here we show a derivation of the formula based on Proposition 15.161 
Note that Mi{p) is common to all the states in one connected component. Hence it is 
determined by the action variable /i and the internal symmetry 7. 

To avoid double indices we denote by f ^ the Zj-th column of the matrix F f l5.32p . 
hence F = {f-^, . . . , fg). For any b & 1^^ we define 

F46] = (/i,...,/,_i, 6, (5.34) 

Let p G VL,-y{fJ^) be a state of the periodic BBS, and S(p) ~ Z^/F-^Z^. We define the 
least common multiple of nonzero rational numbers ri, . . . , r„ as 

LCM(ri,...,r„) = min |Z n Zri n ■ • ■ H Zr„ \ {0}|. (5.35) 



Integrable structure of box-ball systems 57 

For instance we let LCM(l/3,2/3) be 2 rather than 2/3. Now we have 
Proposition 5.18 [571 Th.4.9] 

M(p)^LCm( 't^l.l — h'T^I.i )- (5.36) 

V 7n det Fi [hi] det Fg [hi] J 

where we exclude any entries of LCM such that det Fi[hi] = 0. 

Proof. Proposition 15.161 imphes that N'i{p) is defined as the smallest positive integer 
such that the following linear equation has an integer solution n = (ni, . . . ,ng) G Z^: 

Nhi = F^n. (5.37) 

By demanding all these conditions on the expression Uj = N'ji. det Fj[hi]/ det F of the 
solution to ( ]5.37p . one obtains the formula ( ]5.36p . □ 



Example 5.19 Take p = 1212111222,p = 1211121222 G Vio m Example \5Al\ Then 



det F = det ( ^ 6 ' " "'^^ 



detFi[Aii] = det = 4, detFaf/ii] = det = 4, 

detFi[/i2] = det = 2, detFaf/is] = det = 

det Fi[/i3] = det r g ) = 0, det Fai^s] = det I ^ 3 ) = ^^'^^^ 



Hence 



Ui{p) = LCM(40/4, 40/4) = 10, ^(p) = LCM(40/2, 40/12) = 20, 

7V3(p) = LCM(40/20) = 2. (5.39) 

For p one has 71 = 2,73 = 1. Hence 

7Vi(p) = LCM(40/8,40/4) = 10, N'2{p) = LCM(40/4, 40/12) = 10, 

Afsip) = LCM(40/20) = 2. (5.40) 

For instance, one can deduce T2^{p) = p and T2^{p) 7^ p. One can easily check these 
results by using the figures in ( (5. 33\) . 

When the internal symmetry is trivial, i.e. V7j^. = 1, Proposition 15.181 reduces to the 
result in [831. 
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5.3.5. Multiplicity of torus and structure of isolevel set. The number of the elements 
of the isolevel set Vl^^J') is expressed by the formula [5] 



P.(/^)l = detFn-f^^ + "7') (5.41) 



nr • (5,42) 

Pi \ rrii J 

where the g'xg matrix F is defined as (15.321) . While the former expression (I5.4ip was first 
obtained in the context of Bethe ansatz |^ Th.3.5], the latter one (I5.42p was originally 
found as a formula for the cardinality \Vl{ij)\ and proved by elementary combinatorial 
arguments [831 Prop. 2. 2]. Their equivalence is due to the relation det F = pi^- ■ ■pi^_^L. 

We demonstrate the decomposition of the isolevel set Viif^) into connected 
components from the viewpoint of their cardinality. First we consider the case rrii = 1 
for all i eX. Then by (15.411) one has |'PL(/i)| = detF. Actually this is a consequence 
of Vl{ij) — Jhilj) and (15.191) . In this case Vl{ij) itself is a connected graph. Next we 
consider general cases with rrij > 1. Let m,p be a pair of positive integers and 7 any 
common divisor of m and p. We define 



/ p+m _ 1 \ 

C,{m,p) = Y^J^)[ ^_ , (5.43) 




where /3 runs over all the common divisors of m and p that is a multiple of 7. Here /i 
is the Mobius function in number theory defined by 

if A; = 1, 

/^(^) = \ {^'^y if k is the product of j distinct primes, (5.44) 
otherwise. 

(This yU should not be confused with the Young diagram.) By the Mobius inversion 
formula we have 

7 

where 7 runs over all the common divisors of m and p. By (I5.4ip and (15.450 we obtain 

\VL{^^)\ = ^ mult (7) detF^, (5.46) 
7 

mult(7) = TT^^^^^^^^' (5.47) 

where the numbers {'~fi/mi)C^.{^i^Pi) turn out to be integers as a result of certain cyclic 
group actions [71]. This formula is a consequence of the following fact. 
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Proposition 5.20 [71] The mult(7) in ( [5.^7p is the multiplicity of the tori in Vl{^)- 
That is, the following relation is satisfied: 

~ |Jmult(7) (ZVF^Z^) . 



Example 5.21 Take 7^24 (yw) with 



(5.48) 



Then we have (mi, m2, rn-s) = (3, 2, 1) and (pi,P2,P3) = (12, 6, 4). The matrix F is given 
by 

/l8 4 2 \ 



V 



6 14 4 . (5.49) 

6 8 lOy 

Possible internal symmetries are 7 = (1, 1, 1), (1, 2, 1), (3, 1, 1), and (3,2,1). By using 
^5.43^ and ^5.4T\ ) one has 



P24(/i) ^ 90 (ZVF(i,i,i)Z3) U 30 (ZVF(i,2,i)Z') U 3 (ZVF(3,i,i)Z3) U (ZVF(3,2,i)Z')- 

(5.50) 

The configuration fi is common to that in Example \5.15[ Hence we can deduce that 
S(p) ~ (ZVF(3,i_i)Z3) forp= 121122111212211222121111. 



6. Approach by tropical geometry 
6.1. Preliminary of tropical geometry 

Tropical geometry is the algebraic geometry of the min-plus algebra (T, min, +). We 
introduce the basic notion of the theory of tropical curves, following 



6.1.1. Tropical curve. A tropical polynomial F{X) of one variable X is written as 
F{X) = min(n,X + C,) rii G Z>o, C, G M, 

where / is a finite subset of Z. One can regard F(X) as the tropicalization of a 
polynomial f{x) = ^j^j CiX^' G M>o[2;] with the transformation Cj = e~^ and x = e~T. 
In the same manner, a tropical polynomial F{X, Y) of two variables X and Y is written 

as 

F(X, Y) = min(niX + m^F + Q) Ui, G Z>o, Q G M. 

Very roughly speaking, a tropical curve is a finite graph (i.e. a graph with a finite 
number of vertices and edges) with a metric structure. In the following we only consider 
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Figure 1. Tropical curves: (i) a line, (ii) an elliptic curve 



affine tropical curves in given by a tropical polynomial of two variables. Fix a finite 
set / and a tropical polynomial F{X, Y) of two variables. The tropical curve F given 
by F{X, Y) is defined as 

F = {{X,Y) G I F{X,Y) is mdifferentiable} . 

Example 6.1 See Figure \^ for examples of tropical curves, where (i) is a "tropical 
line" given by F{X^Y) = min(X, Y, 1) and (ii) is a "tropical elliptic curve" given by 
F) = min(2y + 7, 2Y + X + A, Y + 2X, Y + X + 2, F + 6, X + 3, 8). 
The meaning of "indifferentiable" is seen at (i) as follows: let Ai, A2 and 
be three open domains divided by the tropical line. We write I12, I23 and /13 for 
the three boundaries, and P for their intersection point I12 H ^23 H /13. The function 
F{X, Y) = min(X, Y, 1) is "differentiable" at (X, Y) e AiU A2U A3, namely we have 
F{X, Y) = lzn Ai, F{X, Y) =Y m A2 and F{X, Y) = X zn A3. The function F{X, Y) 
is "indifferentiable" at (X, F) G I12 U /23 U /13 where at least two of X , Y and 1 become 
the minimum. For instance, F{X,Y) = Y = 1 on /12 \ {P}, F{X,Y) = X = Y on 
I23 \ {P}, and F{X, Y) = X = Y = 1 at P. 

The edges in tropical curves have rational slopes, and we associate each vertex with 
a primitive tangent vector which is a tangent vector given by a pair of coprime integers. 
(If one of the integers is zero, then let another be ±1.) The primitive tangent vector is 
uniquely determined up to sign. For two vectors ^ = (ni,n2) and = (71^,77,2), we set 
^ A ^' = nin2 — n2n[. The following is the definitions of a smooth tropical curve and its 
genus. 

Definition 6.2 The tropical curve T G is smooth if the following two conditions 
hold: 

(a) all vertices in F are 3-valent. 

(b) For each 3-valent vertex v, let ^1,^2,^3 be the primitive tangent vectors which are 
outgoing from v. Then these vectors satisfy ^1 + ^2 + ^3 = (0, 0) and A = 1 
for i,j G {1,2,3}, i^ j. 
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When a tropical curve T is smooth, the genus ofT is dim ifi(r,Z). 

A smooth tropical curve is equipped with the metric structure as follows: 

Definition 6.3 Assume T is a smooth tropical curve. Let E{r) be the set of edges in 
r, and let be the primitive tangent vector of e ^ E{r). We define the length of edges 
I : E{T) M>o by 

e t-^ Z(e) = -1^4, 

II Se II 

where \\ \\ is any norm in M^. 

Example 6.4 Both of the two tropical curves at Figure U\ 0'1"g smooth. For instance, 
in (a), the 3-valent vertex A has outgoing primitive tangent vectors (1,0), (0,1) and 
(—1,-1), and the conditions (a) and (b) in Definition \ 6.S\ are satisfied. The genera 
of the curves (i) and (ii) are respectively and 1. The lengths of the edges are 
1{AB) = 1{BC) = 1 in (ii), for example. 

We omit the metric structure for non-smooth tropical curves for simplicity. See 



6.1.2. Abelian integral and tropical Jacobian variety. Let F be a smooth tropical curve 
whose genus g is not zero. We fix g generators ■ ■ • , i?^ of the fundamental group of 
r (i.e. generators of the cycles in F). 

To describe the abelian integral on F, we need some preparations: for each e G E{r), 
we fix a linear map : e — )■ [0, 1] (where we have only two choices). For pi,p2 G e, we 
define a fundamental path p by p = (e;pi,p2) ^ E{T) x F x F. For a fundamental path 
p = (e;pi,p2), we define [pi,P2; "e] C [0, 1] by 

[ae{pi), «e(p2)] if ae(Pl) < ae(P2) 
[ae{p2), ae{pi)] if ae(Pl) > 0!e{p2) 

For two fundamental paths p = (e;pi,p2) and p' = {e';p[,P2), we define the addition 
rule only when e = e' and p2 = p[ or pi = p'2, by 

{e;Pup'2) ifp2=p'i 
{e;p[,p2) ifpi=P2. 
We define a set of paths P on F by 

V = [ © Zp) / the addition rule. 

p: fundamental path 

Then V is an infinite dimensional vector space. For two fundamental paths p = (e; pi, P2) 
and p' = {e']p[,P2), we define 

'0 ife^e' 
sgn[{ae{pi) - ae{p2)){ae{p'i) - ae{P2))] if e = e' ' 



p + p 



sgn(p,p') 
lipDp') 



(6.1) 

if e 7^ e' 

\[pi,P2;ae] n [p'i,P2,ae] \ ■ /(e) if e = e' ' 
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and define a bilinear form of fundamental paths by 

( , ) : {p,p') ^ {p,p') = sgn(p,p') - lip Hp'). (6.2) 

Tliis naturally gives the bilinear form ( , ) : V x V ^ M. on V. Briefly speaking, the 
bilinear form of two paths p and p' on F is "the length of the common part of p and p' 
with the sign depending on the directions of the two paths" . 

Example 6.5 See Figure\^for the smooth tropical curve T given by 

F(X, F) = min(2F, F + 3X, Y + 2X, F + X + l, F + 4,11). 

The genus of T is 2, and we fix the basis Bi and B2 of the fundamental group of T as 
depicted. The bilinear forms of Bi and B2 take the values as 

(fii, 5i) = 20, B2) = -7, {B2, B2) = 14. 

Let us demonstrate how to compute {Bi, B2): the common part of Bi and B2 is the edge 
PQ. We set QP C Bi and PQ C B2 which are fundamental paths on T. Then we have 
1{QP n PQ) = 1{PQ) = 7 and sgn(QP, PQ) = -1, and the result follows. 

Now we introduce the abelian integral and the tropical Jacobian variety for F: 
Definition 6.6 Fix Pq ^ T- The abelian integral ip : T is given by 

P^^/;(P) = ((E„P>)),=i,...,3, 

where we choose a path PqP from Pq to P. The map ip induces the map from a set of 
divisors Div(F) on F to M^; 

^ riiPi ^^Ui ^jj{Pi), 
where I is a finite set and rii G Z. 
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Definition 6.7 The period matrix SofTis given by 

E = m,B,)),j=^_,eM,{R). (6.3) 

The tropical Jacobian variety J(T) ofT is the g-dimensional real torus given by 

J{T) = W/EZ3. (6.4) 

Remark 6.8 The matrix S is symmetric and positive definite by definition, and J(T) 
is a tropical analogue of Jacobian variety. By removing all infinite edges ofT, we obtain 
the maximal compact subgraph F'^p* of T. Though the map depends on a choice of 

the path PqP, the induced map T^^^ — j- J{T) does not depend on the choice and becomes 
injective. When g = 1, ip induces the isomorphism from T^^^ to J{T). 

Example 6.9 The tropical curve of genus 1 depicted at Figure U\ (ii) has the period 
matrix S = 9, and the Jacobian variety M/9Z. ^45 for the tropical curve of genus 2 
depicted at Figure\^ the period matrix and the Jacobian variety are respectively obtained 
as 



20 -7 
-7 14 



J(r) = R'/EZ\ 



6.1.3. Tropical Riemann theta function. Fix a positive integer g and a symmetric and 
positive definite matrix S e Mg(M). (Here the matrix S is not always a period matrix 
of some tropical curve.) 

Definition 6.10 The tropical Riemann theta function 0(Z;S) o/ Z G is defined by 

e(Z;S) = minjn- (iSn+Z)}. 

We call the g-dimensional real torus given by 

Jh = MVSZ^ (6.5) 

the principally polarized tropical abelian variety. (IfE. is the period matrix of a tropical 
curve T, then J^ is nothing but the tropical Jacobian variety J(T).) 

It is easy to see the following: 
Lemma 6.11 The function 0(Z) = 0(Z;S) satisfies the quasi-periodicity: 

e(Z + Hm) = -m- (^^Hm + zj +0(Z) m G Z^^. (6.6) 
Remark 6.12 Recall the Riemann's theta function: 

e{z; K) = exp(vri n ■ (A'n + 2z)) z G (6.7) 

where K G Mg(C) is symmetric and \m.K is positive definite. The function 6{z;K) 
satisfies the periodicity and quasi-periodicity: 

e{z + Km- K) = exp(-7ri m ■ {Km + 2z))^(z; K), 
for m G Z^. Remark that only the quasi-periodicity remains in the tropical case. 
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6.2. General solution for tropical periodic Toda lattice 

We briefly present the results on the general solution for the tropicalization of the 
A^-periodic Toda lattice (trop-pToda) following [29]. This section offers not only an 
application of the tropical geometry, but also a preparation for the next section where 
we study the periodic BBS with the similar description to §4.2.21 We emphasize that 
the independent variables of the trop-p Toda in this section are in M, whereas in the 
next section their integer parts (i.e. ultradiscretization) will be studied. 
The trop-pToda is given by the piecewise-linear evolution equation: 



1=1 



(6.9) 



N N 



on the phase space T = {{Qj,Wj)j^z/NZ \ ^Qj <'^Wj} CR^^ . (In [29], this 

system is called the ultradiscrete periodic Toda lattice, where "ultradiscrete" means 
"tropical " in our present terminology, hence the naming "trop p-Toda" here.) 

Proposition 6.13 |42j, Prop. 2.1] Eq. ( \6.9\} is obtained as the tropicalization of the 
discrete Toda lattice (14.181) and (14.191) with a periodic boundary condition (fjj^N ~ 
Qj, ui^j+N — ""^j ^"-^ '^^^ condition Y[f=i'^y < ^ so that the tropicalization of 

Proof. First we show that under A^-periodic boundary condition, (14.181) and (14.191) are 
expressed as 

With a common denominator, the RHS of (16.101) is rewritten as 



Ef=o'ntiK-f+i/gj-;+i) 



On the other hand, (I4.25P are equivalent to (I4.18P and (I4.19P under the substitution of 
variables = l/xj, u;* = l/vj, q^j^'^ = 1/%, w^j^^ = l/vj- Hence the birational R given 
by (14.231) is also equivalent to the A^(= n + l)-periodic discrete Toda lattice equation. 
From (I4.23P with the above variable substitution we obtain 

,+1^ , Ef=inr=.(i/gWi)nti(i/^j-+^) 
''^'EliULiy^umlim+^-i) 
xti uZi'a/QUi^i) nt;(i/^j+.+i) 



Ek=iUi=i\i/QUUi=oiy^u 



^t ^k=illi=k (Wj+i+i/qj^i+i) 
" Ef=inJIfe^K+z/?j+i) 
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which becomes fl6.1ip by the replacements k^N — k,l-^N — I. 

Next we apply tropicalization to flG.lOp and f l4.19p . Since the numerator of the 
second term in f l6.10p is a constant with respect to t, the tropicalization of this numerator 
is constantly zero under the given condition on Y[j'=i ""^j/^j- Thus fl6.10p and fl4.19p yield 

( Km . □ 

The system (16.90 has + 1 conserved quantities {k = 1, . . . , + 1) given by 
the tropical polynomials on T as 

Hi = min {Qj,Wj), 



H2 = mm( min (Qi + Qj), min (Wt + Wj), min (Qi + Wj)), 

l<i<j<N l<i<j<N l<i,j<N,j^i,i-l ' 

N N N 

j=i j=i j=i 



(6.13) 



Fix C = {Ck)k=i,--- ,N+i £ and define the isolevel set Tc by 

Tc = {rer\ Hu{r) = {k = I,- ■ ■ ,N + 1)}, (6.14) 

and the affine tropical curve Tq given by the indifferentiable points of 

F(X, Y) = min(2r, Y + min(iVX, {N - 1)X + Ci, . . . , X + Cn-i, Cn), Cn+i) ■ (6.15) 

We call Tc the spectral curve of the trop-pToda. For the derivation of the conserved 
quantities fl6A3ll and F{X,Y) fl6J[5l) . see [29l §3.1, §3.2]. We set L and A^, rjk for 
A; = 0, . . . , iV - 1 as 

L = Cjv+i-2(iV-l)Ci, 

Ao = 0, Afc = Ck+i-Ck = 1, . . . , iV - 1, 

(6.16) 

rio = L, r/fc = L - 2 ^ min(Afc, Aj) k = l,...,N -1. 

i=i 

The curve Tq is smooth if and only if Ai < A2 < ■ ■ ■ < AAr_i and 77^ > for 
= 1, . . . , — 1. The genus of the smooth Fc is — 1. See Figure 3 for Fc, where we 
set Ci = for simplicity. 

Assume F^ is smooth, and write g for the genus, g = N—1. Fix the basis Bi, . . . ,Bg 
of the fundamental group 7ri(Fc) as Figure 3. In what follows we denote by Q the period 
matrix (16.30 of F^-. The matrix Q = {^ij)ij=i^...^g is explicitly written as 

'?7i_i + r/i + 2(Ai - Ai_i) i = j 

' (6.17) 

-r]j i = J + 1 

otherwise 



and we get the tropical Jacobian variety of Tc as J{Tc) = ]R^/f2Zi^. 

The general solution for the trop-pToda is obtained by the following theorem: 
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Ag-l \g 



Ai A2 

Figure 3. Spectral curve for the trop-pToda 



(6.18) 



Theorem 6.14 When Yc is smooth, we have the following: 

(i) [301 Th. 3.5] Fix Zq G M» and define = Q{Zo + Xt - Leiu), where 

X = (Ai, A2 - Ai, . . . , - A,_i), ei = (1, 0, . . . , 0) G 

The general solution for the trop-pToda is given by 

/^t rpt I rpt-\-\ rpt-\-\ rpt i 

Vn ~ -'n-l "T -'n ~ -'n-l ~ n "r '-^li 

(ii) [nH Th. 1.3] This solution induces the isomorphism JiTc) ~^ Tc- In particular, the 
time evolution of the trop-pToda is linearized on J{Tc), whose velocity is X. 

Example 6.15 The case of N = 2. In this simplest case we can explicitly construct the 
isomorphism a: 

Tc A r^J* ^ JiVc) 

The solution f l6.18p induces the inverse map of ip o a. Let us consider the case of 
C = (0,3,8), where Tc is depicted as 




The following is an example of linearization, where one sees A = (3). We set Pq = O: 
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rc = {{Qi,W,,Q2,W2)} 

t=0 (3,4,0,1) 

t=l (3,1,0,4) 

t=2 (1,0,2,5) 

t=3 (0,2,3,3) 

t=4 (0,5,3,0) 

For general N > 2, the isomorphism 7c J{^c) is regarded as a composition of 
the injective map a : Tc ^ Div^g{Tc) and the abehan integral ip, but a becomes too 
complicated. 

Example 6.16 The case of N = 3 and C = (0,2,6, 19). We have Vt = 
and A = (2, 2). Observe that the velocity in J{Tc) is indeed X. 

rc = {{Qi,W,,Q2,W2,Q3,W3)} ~ J{Tc) 



t=0 (2,1,0,9,4,3) (29,-3) 

t=l (1,0,2,11,3,2) (31,-1) 

t=2 (0,2,4,10,2,1) (10,-10)^(33,1) 

t=3 (1,5,4,8,1,0) (12,-8) 

t=4 (2,7,4,5,0,1) (14,-6) 

t=5 (2,9,4,1,0,3) (16,-4) 



Remark 6.17 Theorem \6.14\ corresponds to a tropical version of J3, {3^ where the 
general solution for the periodic Toda lattice is studied by using (complex) algebraic 
geometry. When Tc is not smooth, neither the structure ofTc nor the solution has been 
clarified yet. It requires a further study on a degeneration of the period matrix Vt and 
Jacobian variety J(Tc)- 

Remark 6.18 In 13^ . AjS and Q were derived by directly ultradiscretizing abelian 
integrals on the spectral curve of the periodic discrete Toda lattice. By combining this 
strategy and tropical curve theory, Theorem \6.14^ i) was proved in IW^ . 

6.3. Periodic BBS and tropical geometry 

As a periodic version of §4.2.21 we have an embedding of the states of periodic BBS 
in those of the trop-pToda ^21 129] . Differently from the case of the original BBS, this 
embedding is not always consistent with the time evolution of the trop-pToda. We 
revisit the results in §5.2.41 with this embedding and tropical geometry. 

Let 'PL(/i) ( 15. 5p be the isolevel set of the L-periodic BBS with the configuration 
11 = {i^'^ , • • • , i^'^)- (Recall we assume zi < ■ ■ ■ < ig.) Set = 1 + Ylk=i ""^h §6.21 



4 4 J{Tc) ^ M/16Z 

(0,7) 9 
(0,4) 12 
(0,1) 15 
(2,2) 2 = 18 

(3, 5) 5 = 21 
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and fix C = (Ci, . . . , Ctv+i) at fl6.14l) as 
Ci = 0, Cn+1 = L, 

_V^- ■ _ _ ^^-^^^ 

Cmi-^+---+mij,_^+«+l — ij'^T^ij + hi ^ — 1, . . . , g , I — 1, 



L, . . . , rriif^. 



Then the embedding 77 : Vl{h) — ^ Tc is defined as follows: among L boxes, fix "the 
leftmost box" (it can be any box) and do the following procedure: 

(i) if the leftmost box is occupied, then set Qi =(the number of the first consecutive 
balls from the left), otherwise set Qi = 0. 

(ii) Set Wi =(the number of i-th consecutive empty boxes from the left) for i = 
I,-- - ,N. U Qi ^ 0, set Qi =(the number of the z-th consecutive balls from 
the left), otherwise set Qi =(the number of the {i — l)-th consecutive balls from 
the left) for z = 2, ■ ■ ■ , iV. 

Then we obtain 

Wi Q2 W2 Wn-1 Qn Wn 

I^^^Il^ 52^T^ when Qi = 0, 

Qi Wi Q2 W2 Wn-1 Qn 

5^1l^£!55T^ when Qi > 0. 

Note that we have Qi = or = 0. 

Example 6.19 The case of L = 9, fi = (3,1) and N = 3. We have C = (0,1,4,9). 
The evolution of the periodic BBS in Viif^i') is at the left, and its embedding in 
Tc = {(Q17 W^i) W^2) W^s)} is in the center. The evolution of the trop-pToda in 
Tc is written at the right. 

t=0 122211211 (0,1,3,2,1,2) (0,1,3,2,1,2) 

t=l 111122122 (0,4,2,1,2,0) (0,4,2,1,2,0) 

t=2 222111211 (3,3,1,2,0,0) (3,3,1,2,0,0) 

t=3 111222121 (0,3,3,1,1,1) (3,1,1,1,0,3) 

The time evolution does not agree with the embedding at t = 3. 

Let Too and TToda be the time evolution operators in Vl{ij) and Tc respectively. As 
one observes in the above example, in general the following diagram is not commutative: 





Tc 






VM A 


Tc 



i.e. ?7 o Too 7^ Troda o V- 
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Proposition 6.20 P^, Prop. 4.4] Let s be a shift operator on Tc given by 
s : {Qi, W,, Q2, W2,..., Qn, Wn) ^ {Q2, W2,..., Qn, Wm, Qi, W^). 

We write Tc/ for the quotient space of Tc with respect to the action of s. Then the 
following diagram is commutative: 

VM ^ Tc/ ~s 
It^ It,... (6.20) 

VM 4 Tc/ ~. 

The induced map rj* gives one-to-one correspondence between ViifJ') and (TcHZ^^)/ 

In Example 16.191 we have s : (3,1,1,1,0,3) t-)- (0,3,3,1,1,1) at t = 3, which 
indicates the commutativity of the diagram fl6.20p . 

Now we come to the final stage of this section. Set rrii^ = 1 {k = 1, . . . ,g), thus 
fl6.19p gives f l6.16p with = ik and rjk = pi^ (15. 8p . Then the spectral curve Tc is smooth, 
and the isolevel set Tc is isomorphic with J{Tc) (Theorem 16. 14p . On the other hand, we 
have the important result on the periodic BBS, the one-to-one correspondence between 
Vl{ij) and J^iif^) due to the map $ (Theorems 15. 8[ ISTTOj) . These results and Proposition 
I6.20l are summarized as follows, which gives a tropical geometrical explanation for J^iifJ')'- 

Proposition 6.21 (^iy)[29l Lemma 2.5] We keep the setting of J{Vc) o^nd Jl{,Ij)- Let 
Jp = W / FIj^ be the principally polarized tropical abelian variety where the matrix F is 
defined by (15.190 . Let c be a shift operator on J{Tc) given by 

c:[(Zi,Z2,...,Z,)]h^[(Zi + L,Z2,...,Z,)]. 

Then we have Jp — J{Tc)/ ~c, where J{Tc)/ ~c is the quotient space of J(Tc) with 
respect to c. 

(a) [291 Eq. (4.5)] We have the following commutative diagram: 

VM ^ Tc/ Tc 

"I" 4-iso. 

JL{^l) c j^^j(rc)/~c ^ J{Tc) 

Example 6.22 We illustrate Proposition 1 6. 21\ (i) by using Example \6.19i The period 
matrices Q (I6.17P and F (15.190 of J{Tc) and Jp are respectively 

Q = 

By changing the basis of iTiiTc) from i?i,i?2 to Bi, Bi + B2, the period matrix Q is 
transformed into 

n' = 

(Of course we have Jn — Jw-) The shift operator c acts on Jqi as a : [{Zi, Z2)] ^ 
[{Zi + 9, Z2 + 9)], and we have R'^/FZ^ ~ J^'/ ~c- TMs indicates Jp ~ JiTc)/ ~c- 
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